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Abstract 

We compute the limiting distributions of the lengths of the longest monotone subsequences of random 
(signed) involutions with or without conditions on the number of fixed points (and negated points) as the 
sizes of the involutions tend to infinity. The resulting distributions are, depending on the number of fixed 
points, (f ) the Tracy- Widom distributions for the largest eigenvalues of random GOE, GUE, GSE matrices, 
(2) the normal distribution, or (3) new classes of distributions which interpolate between pairs of the Tracy- 
Widom distributions. We also consider the second rows of the corresponding Young diagrams. In each 
case the convergence of moments is also shown. The proof is based on the algebraic work of the authors in 
which establishes a connection between the statistics of random involutions and a family of orthogonal 
polynomials, and an asymptotic analysis of the orthogonal polynomials which is obtained by extending the 
Riemann-Hilbert analysis for the orthogonal polynomials by Deift, Johansson and the first author in H. 



1 Introduction 

/3-Plancherel measure 

In the last few years, it has been observed by many authors that there are certain connections between random 
permutations and/or Young tableaux, and random matrices. One of the earliest clues to this relationship 
appeared in the work of Regev |3l| in 1981. A Young diagram, or equivalently a partition A = (Ai, A2, • • • ) h n 
(Ai > A2 > • • • , Aj ; = n) is an array of n boxes with top and left adjusted as in the first picture of Figure 
H|, which represents the example A = (4, 3, 1) h 8. A standard Young tableau Q is a filling of the diagram A by 
numbers 1, 2, ■ • • , n such that numbers are increasing along each row and along each column. In this case, we 
say that the tableau Q has the shape A. The second picture in Figure |l| is an example of a standard Young 
tableau with shape A = (4, 3, 1). Let d\ denote the number of standard Young tableaux of shape A. A result of 
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Figure 1: Young diagram and standard Young tableau 
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The multiple integral on the right hand side is called the Selberg integral, which can be computed exactly for 
each (3 in terms of the Gamma function. In particular, when (3 — 1,2,4, this integral is the normalization 
constant of the eigenvalue density of a random matrix taken from the Gaussian orthogonal ensemble (GOE), 
Gaussian unitary ensemble (GUE), Gaussian symplectic ensemble (GSE), respectively (see e.g. |3(|). Motivated 
by this result, we define the (3 -Plancherel measure on the set Y n of Young diagrams (or partitions) of size n by 

=^3> AGY„. (1.2) 
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A natural question is the limiting statistics of a random A G Y n under M@ as n — > oo. 

The case when (3 — 2 is quite well studied. In this case, M 2 is called the Plancherel measure which arises 

(k) 

in the representation theory of the symmetric group S n . Denote by L n the random variable Afe under the 
Plancherel measure M 2 , and set L n = L„ . In 1977, the limiting expected shape of A under M 2 is obtained in 
[Q, also independently in |}5| for the so-called Poissonized Plancherel measure. In particular, it is shown that 

E(L„) 
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A central limit theorem for L n is then obtained by 
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= F 2 (x), 



(1.3) 



(1.4) 



where F 2 is the so-called Tracy- Widom distribution function, which is expressed in terms of a solution to the 
Painleve II equation (see Definition || in Section || for the definition) . The connection to random matrix theory 
comes from this function F 2 : in 1994, Tracy and Widom fl45[ proved that under proper centering and scaling 
(which is different from the scaling for L n in ( |l.4| )), the largest eigenvalue of a random matrix taken from 
the Gaussian unitary ensemble (GUE) has the same limiting distribution given by F%. In other words, after 
proper centering and scaling, the first row of a random Young diagram under the Plancherel measure behaves 
statistically for large n like the largest eigenvalue of a random GUE matrix. Then in the same paper 0, it 
was conjectured that L„ ^ of a random A € Y n under M 2 have the same limiting distribution as the k th largest 
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eigenvalue of a random GUE matrix for each k. This conjecture was supported by numerical simulations of 
Odlyzko and the second author, and was proved to be true for the second row 1$ in B. The full conjecture for 



the general row L„ was subsequently proved by |38|] , ||10[ and [|32j , independently. The authors in |3S], |1C|, 32 
proved the convergence in joint distribution for general rows, and also in |h], [32|, the authors obtained discrete 
sine kernel representations for the so-called bulk scaling limit of correlation functions, an analogue of the sine 
kernel that appears in the GUE matrix case. In || and ||, in addition to convergence in distribution, the 
authors also proved convergence of moments for L„ and Ln \ respectively. Convergence of (joint) moments for 
the general rows is obtained recently by B . We also mention that there are many works on similar relationships 
between tableaux/combinatorics and GUE random matrices. See for example, 0, |9[ [H], fl4| |28|, |5J We 
refer the readers to [jl], |37], |l2| for a survey and history of L n , and to fl3(| for general reference on random matrix 
theory (see also the recent book Q ) . 

One of the main topics in this paper is the limiting statistics of A 6 Y n under M„. From (1.1) and the 
results for the case when (3 = 2, one might guess that for (3 = 1 the limiting statistics of a random A is same as 
the limiting statistics of the eigenvalues of GOE (Gaussian orthogonal ensemble) matrices. We establish this 
fact for the first two rows. More precisely, we prove (see Theorems ^4 and |3.6| together with the remark that 
follows) , denoting by the random variable Afc of a random A £ Y n under , 

UmPr / £g^Vrc < x)=Fi k \x), k = l,2, (1.5) 



n l/6 

where F±(x) is the limiting distribution function |Q for the (scaled) k th largest eigenvalue of a random 
matrix taken from GOE. We also prove convergence of moments. As in the case of j3 = 2, we expect that the 
above result should extend to the general rows k > 3, and also to the joint distributions. For general values of 



(3 > 0, again from (1.1), we expect that in the large n limit, the rows of a random Young diagram under M@ 
correspond to the Coulomb charges on the real line with the quadratic potential at the inverse temperature 
(3, which specializes to GOE, GUE, GSE eigenvalue distributions for the cases (3 = 1,2,4, respectively. This 
conjecture seems natural from the perspective of the discrete Coulomb gas interpretation for the Plancherel 
measure case M% given by Johansson |3l], |32| . 



Random involutions 

The Plancherel measure has a nice combinatorial interpretation. The well-known Robinson-Schensted 
correspondence [M establishes a bijection between the permutations tt of size n and the pairs of standard 
Young tableaux (P,Q) where the shape of P and the shape of Q are the same and the shape of P (or Q), 
denoted by A(tt), is a partition of n. Thus the Plancherel measure M^ t on Y n is the push forward of the uniform 
probability measure on S n . Moreover, under this correspondence, Ai(-7r) is equal to the length of the longest 
increasing subsequence of tt. More generally, a theorem of Greene Q says that Ai(tt) + •• • + Afc(7r) is equal 
to the length of the longest so-called fc-increasing subsequence of tt. Thus the difference of the lengths of the 
longest /c-increasing subsequence and the longest (k — l)-increasing subsequence of 7r € S n under the uniform 
probability measure is equal to A& of A £ Y n under the Plancherel measure in the sense of joint distributions. 
Thus for example, (O]) and (|l.4|) can be restated for the results on the longest increasing subsequence of a 
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random permutation. On the other hand, the sum of the lengths of the first k columns of A is equal to the 
length of the longest k- decreasing subsequence of corresponding tt. But since the transpose A* have the same 
statistics as A under M%, the results ([ll]) and ( |l.4| ) also hold for the longest decreasing subsequence of a random 
permutation. 

The measure M„ also has a combinatorial interpretation. If tt is mapped to (P. Q) under the Robinson- 
Schensted correspondence, then 7r _1 is mapped to (Q,P) (see e.g. Section 5.1.4 of p3[). Therefore, the set of 
involutions tt = tt^ 1 G S n is in bijection with the set of standard Young tableaux whose shapes are partitions 
of n. Consequently, the uniform probability measure on the set of involutions 

S n = {ir e S n : n = iT^ 1 } (1.6) 



is pushed forward to the 1-Plancherel measure M 1 on Y n . Thus the result (1.5) for k = 1 implies that in the 
large n limit, the length of the longest increasing (also decreasing) subsequence of a random involution behaves 
statistically like the largest eigenvalue of a random GOE matrix. 

An involution tt G S n consists only of 1 cycles and 2 cycles. It turns out that if we put a condition on the 
number of 1 cycles (or fixed points) of tt, the limiting distribution is different. Introduce a new ensemble 

S n , m = {tt G §2n+m ■ \{x ■ n(x) = x}\ = m}. (1.7) 

For an involution tt, the number of fixed points is equal to the number of odd parts of A* (see p3|]). Equivalently, 
the number of fixed points of tt is equal to Ai — A2 + A3 — • • • . Thus the uniform probability measure on the set 
S n ,m is push-forwarded to the measure 

AeU, (1.8) 



where 



{A = (Ai, A 2 , ■ • ■ ) G Y 2n+m : ^(-1)^-% = to}. (1.9) 



j 
j 
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Note that the rows and columns of A G Y n ,m now have different distributions. We denote by L n 'm and L n 'm 
the random variables given by the lengths of the k th row and the k th column of a random A G Y n . m under the 
measure (1.8), respectively. We also set L^ m = Ln,m' and m — L s '^, the length of the longest increasing 



and decreasing subsequences of a random tt G S n ,m under the uniform probability measure. 
Set 



f2n 



(1.10) 



The limiting distribution of differs depending on a. Indeed we prove in Theorem 3.1 and Theorem 9.2 below 
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F 4 (z), 
erf(cc), 



< a < 1, 



a = 1, 



a > 1, 



(1.11) 
(1.12) 
(1.13) 



where F4. and i 7 ! are the distributions for the limiting fluctuations of the largest eigenvalues of random GSE 
and GOE matrices respectively, and erf is the standard normal distribution. Again we also prove convergence 

(2) 

of moments. We note that F4 = F{ ; the limiting distributions of the largest eigenvalue of GSE and the second 
largest eigenvalue of GOE are the same (see the discussion at the end of Section || below.) 

The role of the number of fixed points for the limiting distribution can be seen from the following point 
selection picture. Consider a unit square [0, 1] x [0, 1] in the plane and set 8 — {(x, x) : < x < 1}, the diagonal. 
Suppose we select n points at random in the lower triangle < x < y < 1 and take the mirror image of the 
points about the diagonal 8. We also select m points at random on the diagonal 8. Hence there is a total of 
2n + m points in the square. As illustrated in Figure 0, one such choice of points gives rise to a permutation it 




1 12 3 4 5 6 7 1 
V 1 5 3 7 2 6 4/ 



Figure 2: Point selection process 

satisfying it 2 = 1 with m fixed points, i.e., ir £ S n , m - The length L^(tt) of the longest increasing subsequence 
of 7T is then equal to the 'length' of the longest (piecewise linear) up/right path in the square from (0,0) to 
(1,1), where the 'length' of a path is defined by the number of points on the path. The length of the longest 
up/right path in the above point selection process has the same distribution as L®. Now note that the points 
on 8 form an increasing path. When m is large comparing to n, there are many points on 8 and we expect that 
the longest path consists mostly of diagonal points. Hence we are in the linear statistics situation, and thus 
the order of fluctuation of the length of the longest path is expected to be (mean) 1 / 2 by the usual central limit 
theorem. On the other hand, when m is small comparing to n, then the longest path contains few diagonal 
points (none if m = 0) and we are in the situation of a 2-dimensional maximization problem. In this case it has 
been believed, and in a few cases (for example, [}2| |25|) it is proved, that the fluctuation has order (mean) 1 / 3 . 
(For random permutations, which have a similar interpretation as a point selection process, one can see from 
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the scaling in ( |l.4| ) that the fluctuation is of order (mean) 1 / 3 .) Thus there must be a transition of the limiting 
distribution as the size of m varies. The results (|TTTl) -(p^) show that a = 1 is the transition point. The fixed 
points play the role of adding a special line in the 2-dimensional maximization problem. (See |(| for a relevant 
work where two special lines are added to a 2-dimensional maximization problem.) We note that when a = 1, 
h/5nl * n (Lii) nas tne same limiting distribution as (see (L5)). This is because the typical number of 



fixed points of a random involution of size k is \fk. Indeed the result (1.5) is proved by using (JlTTT|)-([TTT^) and 
taking a summation over the number of fixed points (see Section |^ below) to which only a — 1 gives the main 
contribution. 

Once the transition point a = 1 is known, it is of interest to investigate the transition more carefully. We 

set 

2w 



a = 1 



(2n)V6' 



(1.14) 



and take n — > oo while keeping w fixed. We prove that (Theorem 3.2 below) there is a one-parameter family 
of distribution functions F a (x;w), which is expressed in terms of the Riemann-Hilbert representation for the 
Painleve II equation (see Definition [| below) , such that 



lim Pr 
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n.[a\/2n] 



2V2T 



< x = F m (x;w). 



(1.15) 



(2n + m)V6 

The new class of distributions F a (x; w) interpolates F4 and F± as w — > 00 and w = respectively, and satisfies 
lim^^-oo F a (x; w) — 0, so ( |l,15| ) is consistent with ( 1.11 )— ( 1.13 ). Alternatively, since F4 = F^ 2 \ F a (x;w) 
interpolates the limiting distributions of the second and the first eigenvalues of a random GOE matrix. 

The meaning of F a (x; w) in terms of random matrices is not clear, but there is a Coulomb gas interpretation. 
In (7.64)-(7.65) of (7], the following density function is introduced. Suppose 2N ordered particles on the positive 
real line, < < ■ ■ ■ < £2 < £ii are distributed according to the density function 

2N 



1 



Z2N 



-0/2 



(1.16) 



l<i<j<2N j = l 

where Z2N is the normalization constant. Hence in addition to the usual Coulomb gas interaction, there is 
a additional attraction between neighboring pairs £,2j-i and ^2j, j = 1) ■ • • ,N. When A = 0, this attraction 
vanishes, and one sees that (1.16) is the eigenvalue density for the 2N x 2N Laguerre orthogonal ensemble 
(LOE). On the other hand, when A — > 00, the neighboring particles £,2j, j = lj 2, • • • , JV, coalesce, and 

(1.16) becomes the eigenvalue density function for the TV x N Laguerre symplectic ensemble (LSE). Thus ( 1.16| ) 
interpolates LOE and LSE eigenvalue distributions. This density function arises in a symmetric version of the 
growth model with a growth rule given by the exponential distribution considered in Proposition 1.4 of |pl[ . A 
discrete version of the above density function was considered in (4.27) of S, and the limiting distribution of 
the largest particle was precisely F a (x; w). Now formally taking the exponential limit (set q = 1 — 1/L and 
a = 1 — A/L and take L — > 00) of (4.27) of || (which is convincing, but is not justified yet), we see that if we 
set 



ATl/3 



(1.17) 
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in (1.16) and take N — ► oo, the scaled largest particle (£1 — AN) / {2N 1 / 3 ) has the limiting distribution F a (x; w). 
We plan to exploit the justification of the exponential limit in a later publication. 

On the other hand, the longest decreasing subsequence corresponds to the longest down/right path from 
(0, 1) to (1, 0) in the above point selection process. Thus it is clear that the distribution of L^ m is insensitive 
to the number of fixed points. See Section]^ for results and discussions. 

The other ensemble we consider is the set of signed involutions. A signed permutation n is a bijection 
from {— n, ■ ■ • , — 1, 1, • • • , n} onto itself which satisfies 7r(x) = — 7r(— x). The limiting distribution for a random 
signed permutation is obtained by p7| ^]. In this paper we consider random signed involutions with/without 
constraint on fixed points and also negated points (tt(x) = —x). See Section || for results. Especially, we obtain 
another one-parameter family of distributions which now interpolates F 2 and F^. Here F^ has the meaning as 
the limiting distribution for the largest 'eigenvalue' of the super-imposition of the eigenvalues of two random 
GOE matrices. We note that in this case, F2 is equal to the second largest 'eigenvalue' of such super-imposition 
(see the discussion at the end of Section ||). 

For convenience of future reference, we summarize various definitions introduced above. By the k th row/column 
of 7r we mean the k th row/column of the corresponding Young diagram under the Robinson-Schensted map. 

Definition 1. Let S n be the symmetric group of n letters and let be the set of the bijections from 



n, ■ ■ ■ 



-2, —1, 1, 2, • • • , n} onto itself satisfying tt(x) = —ir(—x). We define 



o,, 



at 



£H,(fc) ( n 
n,m+,m_ V" 



{ir e S n : 7r = 7T 1 }, 

{it e S 2n +m ■ \{x : 7r(x) = a;} I = m}, 

{neS®:n = 

{tt € Sf n+m++m _ : \{x : n(x) = x}\ = 2m+, \{x : tt(x) = -x}\ = 2m,-}, 
the length of the k th row of tt e S n , L n = L£\ 
the length of the k th row of tt e S®, = Lf- (1) , 

the length of the k row Of TT e s® m , 
the length of the k th column of it e S\ 
the length of the k th row of tt e S% 



- £ 0,(D 



s 

n,m 1 



m+ ,m_ 7 



n,m_(- ,m_ n,m+ ,m_ 



(1.18) 
(1.19) 
(1.20) 
(1.21) 
(1.22) 
(1.23) 
(1.24) 
(1.25) 
(1.26) 



The results of this paper were announced in |^]. Since we completed this paper, there have been two 
applications. One is to random vicious walker models [p3| , ^|, and the other is to polynuclear growth models 
p2l |4l| , Indeed there are bijections between the above two applications and various ensembles considered 
in this paper, and thus the results in this paper can be employed to answer asymptotic questions in the above 
applications. 

The proofs of our theorems use the Poissonization and de-Poissonization scheme of ||] . We define the 
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Poisson generating function, for example, for m by (see Definition ^|) 

QfiMM^e-^ £ ^TMC^)' (1-27) 

ni ,n2>0 

A generalization of the de-Poissonization lemma due to Johansson |3(| yields that Pr(L^ 2 nj < I) ~ Q; (rii, ri2) 
as ni, ri2 — ► oo ; see Section |^ for the precise statement. Thus if we obtain the asymptotics of the generating 
function, the asymptotics of the coefficients will follow. The point of the scheme is that the Poisson generating 
functions can be expressed in terms of Toeplitz and/or Hankel determinants. The necessary algebraic work 
for this purpose was done in our earlier paper |Q. The general theory of orthogonal polynomials then tells 
us that Toeplitz/Hankel determinants can be expressed in terms of orthogonal polynomials. In the paper 0, 
it turned out that for all the ensembles being discussed in this paper, we need only one family of orthogonal 
polynomials : Tr n (z;t) = z n + • • • which is orthogonal with respect to the weight e l cos 9 d9 / (2n) on the unit 
circle. This orthogonal polynomial is precisely the same orthogonal polynomial used in || to analyze the 
random permutation problem. The authors in Q| computed the uniform asymptotics of the normalization 
constant N n (t) of 7r„ as n, t — > oo using the steepest-descent analysis for the corresponding Riemann-Hilbert 



problem (see (5J3) below). The difference between the present paper and |3j is that we need 7r„(— a;t) for all 
a > 0, in contrast to || where only one quantity, N n (t), is needed. But in order to analyze N n (t), || controlled 
in a uniform way the asymptotic behavior of the solution to the associated Riemann-Hilbert problem. Therefore 
the asymptotics of 7r„(— a; t) for a uniformly apart from 1 can be (almost) directly read off from the analysis of 



[|[, which eventually imply (1.11). The point z = — 1 (a = 1) in the complex plane plays a special role in this 
Riemann-Hilbert analysis as discussed in || — it is the point where a gap starts to open up in the support of the 
associated equilibrium measure as the relation of t to n varies. When a — > 1 as n — > oo according to the scaling 



(1.14) (which is required for (1.15)) , we need more careful analysis of the Riemann-Hilbert problem which 
is the new part of the asymptotic analysis of the orthogonal polynomials and the Riemann-Hilbert problem. 
In this paper we establish this goal by extending the analysis of B . In the analysis of the Riemann-Hilbert 



problem in Section 10, we give a rather sketchy presentation for the parts which overlap the work of |3j, but we 
give a full proof for the new analysis required for the case a — > 1. We also rework portions of || as necessary 
for consistency of presentation. 

This paper is organized as follows. Section || defines the Tracy- Widom distribution functions as well as 
new classes of distribution functions which will be used to state the main results, and discusses their properties 
(Lemma 2.1). The main results of the paper are then stated in Section ^| Determinantal formulae and orthogonal 



polynomial expressions for the Poisson generating functions are taken from and summarized in Section |]. In 
Section ^|, we state the main estimates of the relevant quantities of orthogonal polynomials; these estimates are 
key to the proofs of the theorems of Section [|. The de-Poissonization lemmas are stated in Section ^. Proofs 
of the main theorems are given in Section [?] for involutions with a constraint on the number of fixed points 



(Theorems 3J, [^|, |3.3| and 3.5), and in Section ||for general involutions and equivalently for (Theorems 3^ 
and [[!]), respectively. The case when a > 1 is considered in Section ^| (see Remark to Theorem |3.3| ). Finally 
the Riemann-Hilbert analysis is given in Section [Io|, which proves the propositions in Section 
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Notational remarks. The ensemble S n in the present paper is identical to in 0. In |Q, is introduced 
to denote the ensemble of neginvolutions, and the authors investigated the longest increasing subsequence of 
7r 6 S„. But there is a bijection between S„ and and the longest increasing subsequence of 7r G S„ 
corresponds to the longest decreasing subsequence of the image of ir. In the present paper, we choose the view 
point of considering both the increasing and decreasing subsequences of involutions of the same ensemble rather 
than considering only the increasing subsequences of involutions of the different ensembles. 

Acknowledgments. The authors would like to thank Percy Deift for helpful discussions and encouragement, 
especially for his help in proving Lemma 2.1. We would also like to acknowledge many useful conversations and 
communications with Peter Forrester, Kurt Johansson, Charles Newman and Harold Widom. Special thanks is 
due to the referee who gave us crucial advice, improving the exposition of the paper significantly. The work of 
the first author was supported in part by a Sloan Doctoral Dissertation Fellowship during the academic year 
1998-1999 when he was a graduate student at Courant Institute of Mathematical Sciences. 



2 Limiting distribution functions 

Let u(x) be the solution of the Painleve II (PII) equation, 

u xx — 2u 3 + xu, (2-1) 

with the boundary condition 

u(x) ~ — Ai(x) as x — > +oo, (2-2) 

where Ai is the Airy function. The proof of the (global) existence and the uniqueness of the solution was first 
established in |27|| : the asymptotics as x — > ±oo are (see e.g. [^7], |l9| ) 

' e -(4/3)x 3 / 2 

Recall that Ai(x) ~ e 2 ^ Fa;1/4 as x — > +oo. Define 

v(x) := { (u{s)) 2 ds, (2.5) 

J oo 

so that v'(x) — (u(x)) 2 . 

We can now introduce the Tracy- Widom (TW) distributions. (Note that q := —u, which Tracy and Widom 
used in their papers, solves the same differential equation with the boundary condition q(x) ~ +Ai(x) as 
x — * oo.) 



u(x) = — Ai(x)+Oy- J, as a; — > +00, (2.3) 

asu — oo. (2.4) 
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Definition 2 (TW distribution functions). Set 



exp 



F(x) 

E(x) := expQ^ 



v(s)ds 

i 

u(s)ds 



exp 



(s — x)(u(s)) 2 ds 



and set 



F 2 (x) :-- 
Fi(i) := 



^(z) 2 
F(x)E(x) 



exp 



(s — x)(m(s)) (is I , 



(F 2 (x)) 1/2 eHM»)^ 



F 4 (x) := F(x) [Eixy 1 + E{x)] /2 = (F 2 (x)) 



1/2 



/2- 



(2.6) 
(2.7) 

(2.8) 

(2.9) 
(2.10) 



In |4q| and [}l6f , Tracy and Widom proved that under proper centering and scaling, the distribution of the 
largest eigenvalue of a random GUE/GOE/GSE matrix converges to F2(x) / Fi(x) / Fi(x) as the size of the 
matrix becomes large. We note that from the asymptotics (2J3) and (2T), for some positive constant c, 

(2.11) 
(2.12) 
(2.13) 
(2.14) 



from the fact that Fp(x) is the limit of a sequence of distribution functions. Therefore Fpfa) is indeed a 
distribution function. 

Definition 3. Define xgoe, Xgue and xgse to be random variables whose distribution functions are given 
by Fi{x), F2(x) and F^x), respectively. Define Xgoe 2 to be a random variable with the distribution function 
F^xf. 



F(x) 


= l + 0(e- cx3/2 ) 


as x - 


■+ +oo, 


E(x) 


= l + 0(e~ cx3/2 ) 


as x - 


-> +oo, 


F(x) 




as x - 


-> — 00, 


E(x) 


= 0(e-^ 3/2 ) 


as x — 


> — CO. 


(x) = 


1 and lim^-oo Fp(x) 


= 0, 


= 1,2,4 



As indicated in Introduction, we need new classes of distribution functions to describe the transitions from 
Xgse to xgoe and from xgue to Xgoe 2 - First we consider the Riemann-Hilbert problem (RHP) for the Painleve 
II equation |2(], [2£j. Let T be the real line R, oriented from +oo to — oo, and let m(- ; x) be the solution of the 
following RHP: 



771(2; x) 

m + (z;x) = rn_(z;x) 




_ e -2j(f z 3 +rrz) ' 



for z e r, 



(2.15) 



l(z;x) = I + 0(~) as z - 

Here m+(z;x) (resp,. m_) is the limit of m(z';x) as z' — > z from the left (resp., right) of the contour T: 
m±(z;x) = lim e ^ m(z =p ie;x). Relation (2.15) corresponds to the RHP for the PII equation with the special 
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monodromy data p = -q = l,r = (sec @ ||, also g§ In particular if the solution is expanded at 



Z = DO . 



we have 



m(z x) = I H h Oi-^j, 

2i(mi(x))i2 = -2i(mi(a;))2i 
2i(mi(x)) 2 2 = -2i(mi(a;))ii 



as z — > oo, 

= u(x), 
= u(ai), 



(2.16) 

(2.17) 
(2.18) 



where u(x) and i> (a;) are defined in ( 2. 1 )-( p.5| ) . 

Now we define two one-parameter families of distribution functions. 



Definition 4. Let m(z;x) be the solution of RHP (2.15) and denote by rrijk(z;x) the (j/c)-entry of m(z;x). 
For w > 0, define 

F a (x; w) :— F(x)l^ [m 2 2(-jw; x) — mi2(—iw; x)] E(xy 1 + [m 2 2(-iw; x) + mi 2 (— iw; x)] J5(x)|/2, (2.19) 



and for w < 0, define 



F^(a;;u;) := e i" ,3 - 2 ™ J F 1 (x) 



Also define 



x < [— m2i(— iiu; a;) + mu(- iw, a;)].E(a;) — [wi2i(— iiM;a;) +mn(— w;x)]-B(a;) >/2. 



F H (x;w) := m 22 (— ?u>; a;)F 2 (a:), w > 0, 

F H (a;;u;) := -e^-^m^Hw; a;)F 2 (a:), «j < 0. 



(2.20) 



(2.21) 
(2.22) 



First, F a (x;w) and F H (x;u;) are real from Lemma 2.1 (i) below. Note that F a (x;w) and F M (x;w) are 



continuous at w = since at z = 0, the jump condition of the RHP ( 2.15 ) implies (mi2)+(0; x) = — (mil)- (0; a;) 
and (m 2 2)+(0;a;) = — (77121)- (0; a;). In fact, F a (x;w) and F ®(x]w) are entire in w £ C from the RHP ( |2.15| ). 
From (|2.1lD - (|2.14| ) and ( |2.24| )-(^.27D below, we see that 



lim F lZl (x;w),F M (x;w) = 1, lim F (x; w), F (x; w) = 

x — >-+oo 



x — > — 00 



(2.23) 



for any fixed w G M. Also Theorem 3.2 below shows that F (a;w) and F H (a;;w) are limits of distribution 
functions, implying that they are monotone in x. Therefore, F a (x;w) and F®(x;w) are indeed distribution 
functions for each ioGR, 

Definition 5. Define Xw an d Xw to be random variables with distribution functions F a (x;w) and F M (x;w), 
respectively. 

We close this section by summarizing some properties of m(—iw; x) in the following lemma. In particular 
the lemma implies that F a (x;w) interpolates between Fi(x) and F±(x), and F^(x;w) interpolates between 



F 2 (x) and F^x) 2 (see Corollary 2.2) 
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Lemma 2.1. Let 03 = ( _°i ) , <J\ = (J J), and set [a, b] = ab — ba. 
(i). For real w, m(—iw;x) is real, 
(ii). For fixed w £ R, we have 



m(-iw;x) = (J+ (e C2;3/2 )) 



'1 -,T 
>0 



Siu 3 — 2xw < 



1 



m(—iw; x) 



1 



V2V 1 1 



ID > 0; .X — > +OO, 

w < 0, x —>■ +00, 



J "I ) e (-|^ 3 +™)<T3 e (^(-a : ) 3/2 + V / 2 l « 2 (-^) 1/2 ) CT 3 j w > Q ; X 



(Hi). For any x, we have 

y , ■ , y f ■ , { ±{E(X) 2 +E(X)- 2 ) -E{ X f 

lim m(-iw.x) = lim (T-imi—iw.xicri — , , 
—0+ v ' w^o- y ' \±(-E(x) 2 + E{x)- 2 ) E(x) 2 

(iv). For fixed i»el \ {0}, m(— iw] x) solves the differential equati 



ton 



—m — w[m, 03] + u(x)a\m, 
dx 



where u{x) is the solution of the PII equation (2.1), (2.2) 
(v). For fixed x, m(—iw;x) solves 



-^—m = (— Aw 2 + x)[m, (73] — Awu(x)aim — 2 ( u , \ \m. 
aw 



(vi). For any x, we have 



i(z;x) = o\m(—z;x)<j\. 



Corollary 2.2. We have 





F®(x;0) 


= F 1 (x), 


lim 

w — *oo 


F a (x;w) 


= F 4 (x), 


lim 

— > — oc 


F a {x;w) 


= 0, 




F K (x;0) 


= ftW 2 


lim 


F m (x;w) 


= F 2 (x), 


w — *oo 


lim 


F®(x;w) 


= 0. 



w — » — 00 
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Proof. The values at to = follow from (2.2S ). For w — > ±00, note from the RHP (2.15) that we have 
limz^oo m(z; x) = I. □ 

Proof of Lemma 2.1. Let v(z) = v(z;x) denote the jump matrix of the RHP (HjJ). 

Since u(— z) = v(z) for 

z S M, M(z) := m(— z; a;) also solves the same RHP. By the uniqueness of the solution of the RHP (2.15), we 
have 



m[-z;x) = m(z;x), z<gC\R. (2.38) 

Thus, m(—iw;x) is real for !«£l, thus proving (i). 

By the symmetry of the jump matrix, a\v(—z) ai = v(z), we obtain, by a argument similar to (i), 

aim(—z;x)ai—m(z;x), (2.39) 

which is (vi). 

The asymptotics results (ii) as x — > ±00 follow from the calculations in Section 6, pp. 329-333, of |]l9f . 
For the proof of (iv), define a new matrix function 



f{z-x) := m{z;x)e- iea;> , 



r 



(2.40) 



Then /(• ; x) satisfies the jump condition f + (z;x) = f-(z;x)( \~^) for z S K, and f(z;x)e 3 — ► / as z — > 00. 
Since the jump matrix for f(z;x) is independent of x, f'(z;x), the derivative with respect to a;, satisfies 
f',(z;x) — f'_(z;x)(\ ^ ) , and f'e %B(j3 + i0' fa^e 19 ' 7 ' 3 — > as z — » 00. Hence /'/ _1 has no jump across R, and 
satisfies + iO'fa^f^ 1 — ► as 2 — * 00. If we write m(z;x) = J + mi j x ) + 0(l/z 2 ) as z — * 00, we have 

i0' fa^f^ 1 — izo~3+i[mi, o~§] +(3(z~ 1 ) as z — ► 00. Thus f'f~ 1 is entire and as z — > 00, f'f^ 1 ~ — izo"3— i[mi, 03]. 
Therefore by Liouville's theorem, we obtain 



f'(z;x)(f(z;x)) 1 = -iza 3 - i [7711,(73] 



(2.41) 



Recalling that u(x) = 2i(mi(a;))i2 = — 2i(mi(x))2i in ( 2.17 ), we have [mi, (73] = iu(x)cri. Changing / to m 

(2.42) 



from (EHl is 



m(z; x) = iz[m(z; x), 0-3] + u{x)a\m{z; x). 



This is ( fT2S| ) when z = — ?u>. 

The proof of (v) is very similar to that of (iv), and the detail is left to the reader. We only note that in the 
derivation of (v), we need the identity 

d 



dx 



(2.43) 



which can be obtained from ( §1§ ) by setting m(z; x) = I + + + 0(1/ z 3 ) as z — > 00. 

Finally we prove (iii). Note that lim lu _^ ± m(—iw;x) = m±(0;x). From the jump condition at z = 0, we 



have 



m+(0;x) = m_(0;x) J 



(2.44) 
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Letting z — > 0, Imz > 0, in (vi), we have crim+(0; x)a± = m_(0;x), which together with (2.44) implies that 
m+(0; x) = <TiTO+(0; x)a\ (\ "q 1 ) . Thus we have 



m+(0; a;) 



a(a;) b{x) 
a(x)+b{x) —b{x) 



(2.45) 



for some a(x), 6(x). Also the condition detw(z) = 1 for all z € K implies that detm(z;x) = 1 for all z 6 C \ ' 
and hence we have 



b 2 + 2ab+l = 0. 



Now letting z — > 0, Imz < in (2.42), we obtain 

Thus from (2.45) and ( [2.46 ), b' /b = —u, which has the solution 

b(x) = b(y)e- J y u{s)ds . 



From (2.24) with w = + , we have 

b(x) = (to 12 )+(0; x) — ► — 1, as an +oo. 



(2.46) 



(2.47) 



(2.48) 



(2.49) 



Therefore, 6(ar) = -e^°° ™( s ) ds , which is -E{x) 2 from (|/7|). Now ( |2.46| ) gives a (a;) = ^(^(x) 2 +.E(a:)- 2 ), proving 



(2.28) 



□ 



3 Statement of results 

3.1 Involutions with constraint on the number of fixed points 

Recall (Definition [j] in Introduction) the ensembles S n ,mj Snm + m_ °f (signed) involutions with constraint on 
the number of fixed (and negated) points. We scale the random variables: 



X 

An,m 



An,m 



(2n + m) 1 /6 



(2n + m)Ve ' 



n,m+ ,m_ 



2^/4n + 2m + + 2m_ 



X "' m+ ^ m - : 2 2 / 3 (4n + 2m+ + 2m_) 1 /6 

Theorem 3.1. For fixed a and (3, we have: 



(3.1) 
(3.2) 
(3.3) 



lim Pr(y r at- i < x ) 
lim Pr( x a < x) 
lim Pr(y , /=- , < an 



F 4 (:r), 
0, 



< a < 1, 



a > 1. 



(3.4) 
(3.5) 
(3.6) 
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(3.7) 



lim Pr(x^,[^ Q] W] < *) = F 2 (x), < a < l,/3 > 0, 

n 1 ™o Pr (<[V^],[^]^ a; )=^W 2 ' 

„ 1 ™ Pr( X ^ [v/SaM ^ ] < z) = 0, a > l,/3> 0. 

As indicated earlier in the introduction, as a — > 1 at a certain rate, we see smooth transitions. 
Theorem 3.2. For fixed w £ K and (3 >0, we have 

lim Pr(x^ m < x) = F (ir;u;), m = [\/2n - 2iu(2n) 1/3 ], 

n — >oo * 

lim Pr(xf, m+ , m _ < = m+ = [Vn - 2wn 1/3 ], m_ = 



(3.8) 
(3.9) 
(3.10) 



(3.11) 
(3.12) 



From Corollary 2.2, this result is consistent with Theorem 3.1. We also have convergence of moments. 
Theorem 3.3. For any p = 1, 2, 3, • • • , the fallowings hold. For fixed a and (3, 



J™ E ((<[vW P ) = E ((XG SE T), < a < 1, 
^((xI [V ^ ] ) p )=H(xgoeT), 

}™ o H(xI [ ^ ] t)=H(xgo E ) p ), 

„ 1 ™ E («[v^,[^]) P )= E ((^o^) P ), 
Also for fixed and f3 > 0, 

lim HtxlmY ) = E(( X S) P ), m = [V^ - 2w{2n) l l% 

n — too ' 

hm E((x^ m+ , m J p )=E((x^n, m+ = [^-2wn 1 'X m_ = 



0</3, 

< a < l,/3 > 0, 
/3 > 0. 



(3.13) 
(3.14) 
(3.15) 
(3.16) 
(3.17) 

(3.18) 
(3.19) 



Remark. Theorem 3.1 shows that when a > 1 is fixed, we have used incorrect scaling. When properly scaled, 
the resulting limiting distribution is Gaussian. See Section ^ for the statement and the proof. 

The proof of Theorem 3T, [T2, 3J3 are provided in Section f[[ 

In terms of the points selection process, which is a version of (directed site) percolation, mentioned in the 
Introduction, the above results show that the limiting distribution of the longest path depends on the geometry 



of the domain, while the order of fluctuation is the same: (mean) 1 / 3 . From (1.4), the longest up/right path in a 
rectangle < x, y < 1 has F% in the limit, while the longest up/right path in a lower triangle < x < y < 1 has 



i*4, (3.4), in the limit if there is no points on the edge < x = y < 1. If there are points on < x = y < 1, they 
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affect the length of the longest up/right path. On the other hand, the longest down/right path corresponding 
to L^m can be thought as the longest path from the point (0, 1) to the line < x = y < 1. Thus the result 



(3.7) shows that the point-to- line maximizing path has different limiting distribution from the point-to-point 



maximizing path, F% from (E4), though the fluctuation order is identical. One can also state similar results for 



Poisson process (see Theorem \LQ below) and also certain directed site percolation processes considered by |31 
(see ||). This observation came from discussions of one of us (J.B.) with Charles Newman to whom we are 
specially grateful. 

3.2 General involutions 

Now we consider general involutions and signed involutions without any conditions on the number of fixed or 
negated points. 



Theorem 3.4. For any fixed x £ K, we have 



Jim Pr( xn ■= nl/6 v < x ) = F t (x), (3.20) 



Also for any p = 1, 2, 3, 



hm E((xn) p ) =E(( X goe) p ), (3.22) 

n—too 

hm E((xf) p ) =E(( XGOB2 ) P )- (3.23) 

n^oo 

As mentioned in the introduction, this result proves that the first row of a random Young diagram under the 
1-Plancherel measure M\ behaves statistically like the largest eigenvalue of a random GOE matrix as n — ► oo. 
The proof of the theorem is given in Section |^. 

3.3 Second rows 



For the second row, we scale the same way as in (3.1)-( |Q| ), and denote the scaled random variables by Xnffl \ 
Xn,m ) and Xn',m+,m-, respectively. 

Theorem 3.5. Let a, (3 > be fixed. Then 

lira Pr( X a % ) - Ft(x), (3.24) 



n-^-oo 



16 



and for any p — 1, 2, 3, • • • , 

J™ E (^;SrJ p ) = E ((^) p )' ( 3 - 27 ) 

}^ E ^l%i 0] r)= E ((x G s B n (3-28) 

^((Xn'^u^T) = H(xgue7). (3.29) 
This theorem is proved in Section 

As in the first row, these results yield the following theorem on the second rows of general (signed) involutions. 



The proof is very similar to the proof of Theorem 3T and we skip the details. 
Theorem 3.6. For any fixed we have 



Km Pr( xl 2) := — Jf 1 < x ) = F,(x), (3.30) 



„i/e 

Also for any p = 1, 2, 3, • • • , 

lim E((xl 2) ) p ) =E((xgse) p ), (3.32) 

n^oo 

lim E((x^ (2) ) p )=E(( XGf/B r). (3.33) 

n — >oo 

We conclude this section with some remarks on GOE and GSE. If the conjecture given in the introduction 
that the k th row of a random involution behaves in the limit like the k th largest eigenvalue of a random GOE 
matrix were true, the result (3.3C) suggests that the limiting distribution, F± , of the second largest eigenvalue 
of GOE is equal to the limiting distribution, of the largest eigenvalue of GSE. Equivalently, since a GSE 
matrix has double eigenvalues, the second eigenvalues of GOE and GSE are expected to have the same limiting- 
distribution. An indication for this is Theorem 10.6.1 of |36|], which says that the distributions of N alternate 
angles of the eigenvalues of a random 2A^ x 2A^ matrix taken from the circular orthogonal ensemble (COE) 
are identical to those of the N angles of the eigenvalues of a random N x N matrix taken from the circular 
symplectic ensemble (CSE). Indeed for 2A^x 2A^ Laguerre ensembles, we have proved that the joint distributions 
of the second, fourth, sixth, . . . largest eigenvalues of the Laguerre orthogonal ensemble (LOE) and the Laguerre 
symplectic ensemble (LSE) are identical (see Remark 1 to Corollary 7.6 of 0). In particular, since the k th 
largest eigenvalue of a Laguerre ensemble has the same limiting distribution as the corresponding quantity for 
a Gaussian ensemble, the above remark implies that 

F m-x) = F (2k) = F ( 2fe ) ; k = 1)2 ,.... (3.34) 



Thus (3.2C) and (3.3C) imply that the first and the second row of a random involution has the same limiting 
distribution as the first and the second eigenvalue of GOE, respectively. 

Recently in [p4| , the authors proved that the same property holds true for GOE and GSE. They also proved 
that, among many things, the (2k) th 'eigenvalue' of a superimposition of two random GOE matrices has the 
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same distribution as the k eigenvalue of a random GUE matrix. In particular, when k = 1, this implies that 



(F^=F 2 , 



(3.35) 



and hence (3.21) and (3.31) states that the first and the second rows of a random signed involution have the 
same limiting distribution as the first and the second 'eigenvalues' of the superimposition of two random GOE 
matrices, respectively. 



4 Poisson generating functions 

We review the results from JtJ which we will need in the proof of the theorems of Section [| As in |7j , throughout 
the paper, the notation © indicates an arbitrary member of the set {0, S, 0}. 

Definition 6. We define the Poisson generating functions for the distributions introduced above: 

QF(Ai,A 2 ) := e-*-* £ ^^C^')- (4-1) 

ni ,7i2>0 

Qf(A l5 A 2 ) := £ ^- Pr (^ < Z ) ; (4 . 2) 

ni ,7i2>0 

\™1 \™ 2 \"3 

Qf(Ai, A 2 , A 3 ) := e -^"^->3 £ A t A 2 A 3 pfojj. < f ), (4 3) 

^— ' . n ?ii!n 2 !?i3! 

ni,n2,n3>0 

As in JlJ, let /jf mJ (resp., f^mi) ^ >e the number of involutions on ?i numbers with m fixed points with no 
increasing (resp., decreasing) subsequence of length greater than Thus /( 2 „ 2+ „ 1 )„ 1 ( = P r (^n 2 ,ni — ' \^n2,m |i 
etc. Also let /^f m+m _/ be the number of signed involutions on In letters with 2m + fixed points and 2m_ negated 
points with no increasing subsequence of length greater than I : /p n3+ni+n2 ) nin2 ; = p *( L n s ,n lt n 2 ^ l Y\ S n a ,n u n 3 \- 
We also define 

P?(t;a) := c — "* 2 / 2 E E «"7»™i. ( 4 - 4 ) 

0<n " 0<m 

J.TI 

if fe/J) := e* 2 / 2 E E ( 4 - 5 ) 

0<n ' 0<m 

if (f ; a, 0) := e -«*-/«-* E ^ E C-l' ( 4 ' 6 ) 

0<n 0<m+,m_ 

Using I S n m I — ^ 2 ? + .!?2 ! and IS,?™ ™ I = ( 2n + m ++ m -) ! ^ j s eaS y ^ c heck that 

P?(t;a) = Qf(at,t 2 /2) : (4.7) 
P ; s (i;/3) = Qf(/3t,t 2 /2), (4.8) 
if(i;a,/3) = Qf(at,(3t,t 2 ). (4.9) 



It turns out that the P-formulae in (4.4)-(|4.q) are useful for algebraic manipulations (see |7j]), while the Q- 
formulae (4.1)-(4.3) are well adapted to asymptotic analysis. 
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The following results from provide the starting point for our analysis in this paper. For a nonnegative 
integer k, define Wk{z;t) — z k + ■ ■ ■ to be the monic orthogonal polynomial of degree k with respect to the 
weight function exp(t(z + 1 / z))dz / (2ni) on the unit circle. Let the norm of TTk(z;t) be Nk(t) : 

f i 7 

7r„(z;t)7r m (z;t)e^ +1 / z >-— = N n (t)S nm . (4.10) 
s Zmz 

We note that all the coefficients of 7r„(z;i) are real. Define 

<(z;t):=z\(«- 1 ;t). (4.11) 

Then 

Theorem 4.1 (Collorary 4.3 and Collorary 2.7 of @). For a, 13 > 0. 

P 2 7(i;a) = e- at -* 2 / 2 i|[7r^_ 1 (-Q;f)-a7r 2 i-i(-a;t)] J D ; — (t) + [^(-ajt) + mra-x (-a; t)] £>++(*) J, 

(4.12) 

^+1 (*l = e^ Qt -* 2 / 2 ~ ( [tt*, (-a; t) + cw 2 , (-a; t)] e*D+" (t) + [n* 2l (-a; i) - an 2l (-a; *)] e-^- + (* ) 

(4.13) 

P 2 c f +1 (i;/3)= e - t2 / 2 J D++(t), (4.14) 
Pjf +1 (t;a,/?) = e^'V ("«;*) A(t), (4.15) 

where for any real t > 0, 13; (i) and D t (t) are certain Toeplitz and Hankel determinants which in turn can be 
written as 

e-* 2 A(*) = Il^iW" 1 ' ( 4 ' 16 ) 
e"' 2 / 2 ^— (i) = JJj^+sW-^l+TTiy+aCOjt)), (4.17) 
e~ t2 l 2 D+ + {t) = JJj^+aCtJ-^l-Try+aCO;*)), (4.18) 

e-*'™Dt-{t) = n^+i(*) _1 (l-^+i(05*))> ( 4 - 19 ) 
e - t V 2 - t£) - +(t) = J]iv 2 . +1 ( t )-i(i +7r2i+1 (0;t)). (4.20) 



Remark. The absence of /? on the right hand side of ( 4.14 ) is fairly simple to explain. Observe that in the point 



selection model, the longest decreasing subsequence can always be chosen to be symmetric about the diagonal; 
moreover, any decreasing subsequence can contain at most one diagonal point. Thus if the longest decreasing 
subsequence has I points, then removing the diagonal points will result in a longest decreasing subsequence 
with 2 [i /2] points. The independence from /3 is thus special to P ; s for I odd; for I even, we do indeed have 



/3-dependence. But by the monotonicity of P ; in I, we only need (4.14) to compute the limiting distribution; 



in particular, the limiting distribution will not depend on f3. A similar remark applies to (4.15). 
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As a special case: 

Theorem 4.2 (Theorem 2.5 and Collorary 4.3 of [Q]). For I > 0, we have the following formulae: 

Pl +2 (i;0) = e-* 2 / 2 [Di^W+Dt+Q)] /2, (4.21) 
P| +1 (i;0) = e~* 2 / 2 [e*Z?+-(t) + e-*Ur + W] A (4-22) 
P®(t;0) = e- t2 / 2 D++(t), (4.23) 

if (*; 1) = P$(t\ 1) = e- t - t2 / 2 D I - + (*). (4-24) 

P$ + ifr 1) = Pi +1 {t; 1) = e-*'*D++(t), (4.25) 

P|(i;0,0) = e-* 2 AW, (4-26) 

P^ +1 (t;l,/3)=e- t -* 2 D++(t)£> r +(t) ) (4.27) 

Pj? +3 (t; 1,(3) = e-*-* 2 ^^)^^)- (4.28) 

P o (t;O) = er t2 l 2 D^-{t) = e~* </2 . 

For the second row, we define the Poisson generating functions in a similar manner. Then we have 

Theorem 4.3 (Theorem 5.8 and Collorary 5.12 of @). Fora,f3 > 0, 

P^ 2) (t;a)=P?(t;0), (4.29) 
P 2 ^(*;/?)=P 2 ^;0), (4.30) 
P®£\t;a,(3)=PS(t;0,0). (4.31) 

5 Asymptotics of orthogonal polynomials 

Let £ = {z G C : |z| = 1} be the unit circle in the complex plane, oriented counterclockwise. Set 

^(«;t):=e*(' + *~ 1 ). (5.1) 
Let 7r„(z; t) = z" -I be the n-th monic orthogonal polynomial with respect to the measure i/j(z; t)dz / {2iriz) 



on the unit circle. From Theorems 4.1 and 4.2, in order to obtain the asymptotics of the Poisson generating 



functions, we need the asymptotics, as k,t — > oo, of 

N k (t), n k {z;t), (5.2) 

In this section, we summarize the asymptotic results for these quantities. 
Define the 2x2 matrix- valued function of z in C \ £ by 

Y(z;k;t):= [ H ' J ^^T^M,^ 1, k > 1. (5.3) 



20 



Then Y(- ; k;t) solves the following Riemann-Hilbert problem (RHP) (see Lemma 4.1 in p|): 

Y(z; k; t) is analytic in zeC\E, 



Y + (z;k;t) = Y_(z;k;t) 



Y(z;k;t) 



= / + 0(i) 



as 



(5.4) 



>D 



Here the notation fe) (resp., Yl) denotes the limiting value lim z /_» z Y(z'; k) with \z'\ < 1 (resp.. 

Note that k and t play the role of external parameters in the above RHP ; in particular, the term 0(|) does not 
imply a uniform bound in k and t. One can easily show that the solution of the above RHP is unique, hence 
(5.3) is the unique solution of the above RHP. This RHP formulation of orthogonal polynomials on the unit 
circle given in ]3| is an adaptation of a result of Fokas, Its and Kitaev in El]] where they considered orthogonal 
polynomials on the real line. 



From (5.3), the quantities in (|5.2|) are equal to 



= -y 21 (0;M), 

= Yu(z;k;t), 



Y 21 (z;k + l;t)(Y 21 (0;k + l;t)Y 



(5.5) 
(5.6) 
(5.7) 



(For the other entries of Y, one can check directly from ( |5.3| ) that Yi 2 (0; k;t) = Nk(t), ^22(0; k;t) = 7Tfc(0;i).) 
Thus the asymptotic analysis of the RHP ( |5.4| ) would yield the asymptotics of the above quantities, and hence 
eventually the theorems in Section ||[ The asymptotic analysis of the RHP (5.4) was conducted in || with 
special interest on Y 2 i(0; k;t). But as mentioned in the Introduction, || controlled the solution Y(z) to the 
RHP (5.4) in a uniform way. In || and the theorem below, it is natural to distinguish five different regimes of k 
and t. From the analysis of ||, the following results for Y(0; k; t), except for 7Tfc(0; t) in the case (ii) below, can 
be directly read off. For example, (5.34)-(5.35) of || yield the estimates for the case (iii) below when x > 0. 
For the case (ii), we need to improve the L 1 norm bound, (5.23) in ||, of the associated jump matrix. If one 
is interested only on Nk-i(t), the first integral involving w^(s) in the displayed equation before (5.19) of (|| 
vanishes, and hence the bound (5.23) of || was enough. But for 7Tfc(0;t), this integral does not vanish, and we 
need an improved bound. See the discussion in ( 10.43 )-( 10.45 ) below. 

Proposition 5.1 ([§). There exists Mq > such that as k,t —> oo, we have the following asymptotic results 
for Nk-i(t) and 7Tfc(0;t) in each different region of k and t. 

(i)- If0<2t<ak with < a < 1, then 



|iVfc-i(t)~ 



kfc(0;t)| < Ce 



- ck 



(5.8) 



for some constants C, c > 0. 
(ii). If ak <2t <k- Mk 1 / 3 with some M > Mq and < a < 1, then 



(5.9) 
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for some constant C > 0. 



(in). If2t = k— ^TTak 1 / 3 with —M < x < M for some constant M > 0, the 



Nk-i(t)' 



1 



2 l/3 



k l/3 



7r fe (0;i) + (-l) k ^ u (x) 



< 



C 



for some constant C > 0, where u(x) and v[x) are defined in (2.1) and (2.5) respectively 
(iv). Ifk + Mk 1 / 3 <2t<ak with some M > M and a > 1, then 

iTt 



,fc(2t/fc-log(2t/fe)-l)j Vfe _ i m-l _ J 



2t 
2< - fc 



7r fc (0;t)-l 



< 



C 



2* - jfe 



/or some constant C > 0. 
f'y). 7/afc <2t <bk with 1 < a <b, then 
^2t 



k(2t/k~log(2t/k)-l) Nk _^ t yl _ 1 



2/ 



2t — k 



7r fc (0;t)-l 



< 



C 



/or some constant C > 0. 



(5.10) 



(5.11) 



(5.12) 



Now we are interested in 7Tfc(z;t). If z is apart from —1 and is fixed, then similar estimates for y(z;fc;t) 
can be obtained from the analysis of ||. The result below when x > is (almost) direct from the work of ||. 
For the case when x < 0, the analysis of |3| expresses the bound in terms of the so-called (/-function, and we 
need further analysis for this g-function. When z = 0, this ^-function becomes very simple: g(0) = iri. See 



( |10.134D -( [10.139D below. 

Proposition 5.2. For 2t = k — x(k/2) 1 ' s , x fixed, and for each fixed z € C \ S, we have 



lim e tz TTk(z; t) = 0, 

k — >oo 



lim fe _ 



= 1, 



\z\ < 1, 



lim z-*e tz ~n k (z;t) = 1, lim*-,*, *-*e ta Vjjfot) = 0, |z| > 1 



(5.13) 
(5.14) 



Corollary 5.3. For 2t = — x(k/2) 1 / 3 , cc /lied, we ftave for fixed a > 1, 

lim e" at 7r fe (-a;t) = 0, lim e~ Q *7r£(-a; i) = 0. 

k — >oo — *oo 



Proof. Write 



where 



e- Qt 7r fe (-a;f) = a k e t( ~ a+a >aT fe e-* a ir k (-a;t) 

kf(a-2t/k) -k -ta _1 _ / a 

— e J y ' ' ' a e 7Tfc(— a; r), 



/(a; 7) = -^(~a + a 1 )+loga. 



(5.15) 



(5.16) 



(5.17) 
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The function /(a; 1) is strictly decreasing for a > 0, and /(l; 1) = 0. Hence /(a; 1) < for a > 1. Note that 

(5.18) 



f(a; 1 ) = f(a;l) + ^-l(-a + a- 1 ). 



When x < 0, 2t/fc > 1, hence f(a;2t/k) < f(a;l). On the other hand, when x > 0, since 2t/fc — 1 

\ 3/2 
2 2 /3 x( _ a+a -l^ 



-a;/(2 1 /3fc 2 /3) 5 we have f( a . 2t/k) < \f(a; 1) if k > 



\e- at ir k {-a;t)\ < e 



(-a)- fc e-* Q *7r fc (-a;t)| -> 0, 



Therefore (5.14) implies that 



as k — > co. Similar calculations give the desired result for 7T? 



(5.19) 
□ 



When z — > — 1, (which is required for the proof of Theorem 3.2), the estimates for Y(z;k;t) can not be 
directly read off from the result of ||. However, with more detailed estimates, the same procedure as in (3) 
gives us the following restuls. See Subsubsection 10.1.3| ( |10.73[ ) - ( |10.83| ) for the case x > 0, and Subsubsection 
10.2.3| QlO.140| ) - ( |10.150D for the case x < 0. Recall from Section | that m(- , x) solves the RHP for the PII 



equation ( 2.15] ). 

Proposition 5.4. Let 2t = k — x(k/2) 1 ^ 3 where x is a fixed number. Set 

2 4 / 3 w 



a = 1 



We have for w > fixed, 



lim (— 1) e Q 7Tfc(— a;t) = — m\%{— iw; x), 



and for w < fixed, 



lim e _tt X(-a;i) 



lim {-a)- k e- ta K k (-a;t) 

k — >oo 



lim a e 

k — >-oo 



k _ — ia 



7r£(-a;f) 



m 2 2(-iw; x), 

- mii(—iw;x), 
-- -m2i(—iw;x). 



Corollary 5.5. For w < 0, under the same condition as the above proposition, we have 

lim (-l)- k e- ta TT k (-a;t) = m n (-iw; x)e {s/3)n]3 - 2xw , 



lim e * a 7Tfc(— a; t) = 

k — >oo 

Proof. Note that under the stated conditions, we have 



-m 21 {-iw-x)e^^ w3 ' 2xw . 



^(a-i-a)^ _ e (8/3)w 3 -2xw+0(k-^ 3 )_ 



(5.20) 

(5.21) 
(5.22) 

(5.23) 
(5.24) 

(5.25) 
(5.26) 

(5.27) 
□ 



Remark. As noted in Section |2j, it follows from the RHP ( 2.15 ) that (77112)+ (0; x) = -(mu)_(0;i) and 
(m22)+(0;x) = — (m2i)_ (0; x), and hence by the above Corollary, the limits in (5.21), ( [5.22 ) are in fact contin- 
uous across w = 0. 
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For convergence of moments, we need uniform bound of n k (z;t) for \x\ > M for a fixed number M > 0. 
The results ( 5.29 ) and ( 5.30| ) are essentially in the analysis of S, while ( 5.31 )-( 5.34 ) are new estimates. We 
again need to extend the method of 0] to obtain the results below. The proof is provided in Section [tol See 



Subsubsection 10.1.1 and 10.1.2 for the case x > M, and Subsubscction 10.2.1 and 10.2.2 for the case x < —M 



Proposition 5.6. Define x through the relation 

2t 



(5.28) 



2V3fc2/3" 

Then there exists Mq such that the following holds for any fixed M > Mq . Let < b < 1 and < L < 2~ 3 l 2 \fM 
be fixed. Then as k,t — > oo, we have for x > M , 



\e tz n k (z;t)\ < Ce- C ^'\ 
|e"~V fc 7r fc (z;t)-l| < Ce~ c ^'\ 
c (-a;t)l <Ce c l*l, 



\e~ ta Tr k ( 



\e~ ta ~ {-a)- k ir k {-a-t) - l| < Ce"^ 1 



3/2 



\z\<b, 
\z\>b~\ 

a=l- 2 i ' 3 k- 1 / 3 w, 
a=l- 2^ 3 k^ 3 w, 



-L <w < L, 
-L<w<L, 



(5.29) 
(5.30) 
(5.31) 
(5.32) 



and for x < —M, 



V fc (-a;t)l < C, < a < 1, 



{-a)- k TT k {-a;t)\ < C, 



a > 1. 



(5.33) 
(5.34) 



Corollary 5.7. Let a — 1 — 2 4 / 3 wk 1 / 3 and —L < w < L for fixed L > 0. Under the assumption of the above 
proposition, for x < —M, we have 



\e- ta n k (- a] t)\ < Ce^\ 
-\-a)- k ir k (- a] t)\ < Ce^. 



for some positive constants C , c. 
Proof. We have 



\e t( - a - a ^a~ k \ 



-2xw-^w 3 +0(k~ 1 ) 



Above proposition shows that ( 5.35 ) is true for w > 0. For w < 0, write 



e- ta ir k (-a;t) 



{-a) k 



{-a) k ir k (-a;t) 



Now ( |5.35 ) follows from (5.34) and (5.37). The estimate (5.38) is proved similarly. 



(5.35) 
(5.36) 



(5.37) 
(5.38) 



(5.39) 
□ 



The result below is new and is used for the asymptotics of , r— , when a > 1. See Subsection 10.3 for 

■' * n,[V2na] I 1 

the proof. 
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Proposition 5.8. Let a > 1 be fixed. When 

t _ a a(a 2 - 1) 1/2 x 



k a 2 + l (a 2 + I) 3 / 2 Vfc 



— =, x fixed, 



(5.40) 



lim e^-o^TT^-aT 1 ;*) 

k^oo 



'2-k J- 



e 2y dy. 



(5.41) 



6 De-Poissonization lemmas 



In this section, we describe a series of Tauberian type de-Poissonization lemmas, which enable us to extract 
the asymptotics of the coefficient from the knowledge of the asymptotics of its generating function. Lemma 6.1 
below is due to Johansson ]3(| and Lemma 3/2 is taken from Section 8 in || . Lemma |6.3| and Lemma [0| are 



multi- index versions. Lemma 6.1 and 6.3 are enough for both convergence in distribution and the convergence 



of moments, but for convenience of computations, we use Lemma 6.2 and |6.4| for the convergence of moments 
in the subsequent sections. 



For a sequence q = {q n }n>o, we define its Poisson generating function by 



</>(A) 



y^Qn—r- 

t—^ n,\ 



Q<n 



Lemma 6.1. For any fixed real number d > 0, set 



,(d) 



.+ (2Vd+l + l)Vnloi 



vi d) = n - {2Vd+l + 1) log ra. 



(6.1) 



(6.2) 
(6.3) 



Then there are constants C and n a such that for any sequence q = {q n } n >o satisfying (i) q n > q n +i (H) 
< q n < I, for all n>0, 



4>{^) - Cn- d < q n < 4>{ V ^) + Cn- d 



(6.4) 



for all n > Uq. 



Lemma 6.2. For any fixed real number d > 0, there exist constants C and no such that for any sequence 
q = {q n }n>o satisfying (i) and (ii) above, 



q n < C(p(n — dy/n), 
l-«n < C(l-$(n + dy/n)) 

for all n > Uq. 

For multi-indexed sequences, there are similar results. For q = {(3 , ni,n 2 }ni,n 2 >o, define 

X™ 1 X™ 2 

<!>(\ 1 ,\ 2 ) = e-^ Yl fe-Vr- 

^— ' ni\n 2 \ 
From the above two lemmas, we easily obtain the following lemmas. 



(6.5) 
(6.6) 



(6.7) 
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Lemma 6.3. For any real number d > 0, define fin an d vffi as in Lemma 6.1. Then there exist constants C 
and n such that for any q = {q ni ,n 2 } ni ,n 2 >o satisfying (i) q nit „ 2 > q ni +i.n 2 , Qm,n a > ?ni.n 2 +i (&>) < ?ni,n 2 < 
1, for all ni, n 2 > 0, 

0(mS?,/42) °( n i d + n 2 d ) < Qn in2 < + C{rq d + n 2 d ) (6.8) 

for all n\,U2 > Uq. 
Similarly, 

Lemma 6.4. For any fixed real number d > 0, there exist constants C andno such that for any q = {qm,n 2 }ni,n 2 >0 
satisfying two condition in Lemma 6.5, 



q ni n 2 < C(f>(ni — dy/rH,Ti2 — dy/rio), 
1 - q ni n 2 < C(l - (f>(ni - ds/nl, n 2 + dy/ivz)) 

for ni,n 2 > n . 

Remark. Similar lemmas hold true for sequences of arbitrarily many indices. 



(6.9) 
(6.10) 



7 Proofs of Theorems |37T], gg, and ^ 



The following results follow from Proposition 5.1 above. The result ( |7~l| ) is derived in Lemma 7.1 (iii) of 
and the other cases are similar. We omit the details. 



Corollary 7.1. Let M > Mq, where Mq is given in Proposition 5.L. Then there exist positive constants C and 



c which are independent of M, and a positive constant C(M) which may depend on M, such that the following 
results hold for large I. 

(i)Bj Define x by 2t = I — x(l/2) 1 / 3 . For —M < x < M, 



^log^(<)- 1 -21ogF( a 

3>l 



(ii) Define x by t = I - (x/2)l 1 / 3 . For —M < x < M, 

^log^+a^)- 1 -\ogF{x) , ^logjy^+i^)" 1 -logF(x) 

3>l 3>l 

^log(l-7r 2i+2 (0;t))-log£(x) , £ log(l + 7r 2j - +1 (0; t)) - logE(x) 



3>l 



3>l 



^log(l+7r 2i+2 (0;t)) + log£(a;) , ^log(l - 7ny+i(0;t)) + \ogE[x) 

3>l 3>l 



< 



< 



< 



C(M) M 3/2 

£1/3 " 1 "'~ /e ' 

C{M) 



11/3 

cm 

11/3 



■Ce 



Ce 



(7.1) 



(7.2) 
(7.3) 
(7.4) 



These results yield the asymptotics of the determinants in Theorem 4.1 
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Corollary 7.2. There exits M\ such that for M > Mi, there exist positive constants C and c which are 
independent of M, and a positive constant C(M) which may depend on M such that the following results hold 
for large I. 

(i) Define x by 2t = I - x(l/2) 1 / 3 . For -M < x < M, 



le-' Di(t)-F(x) 2 \ < 



C(M) 
1W 



Ce 



cM 3/2 



(ii) Define x by t = I - (x/2)l 1 / 3 . For —M < x < M, 

\e- t2/2 D--(t)- F(x)E( X y 1 \ < 

\e- t2/2 D++(t) ~ F(x)E(x)\ < 

\e- e/2+t D+-{t)- F{x)E(x)- 1 \ < 

\e- e/2 - t Dl+(t)-F{x)E{x)\ < 



C(M) 



11/3 

C(M) 
C(M) 



+ Ce- 



cM 3 ' 2 



Ce 



Ce 



-elf 3 ' 2 



cM 3 ' 2 



(7.5) 

(7.6) 
(7.7) 
(7.8) 
(7.9) 



Proof. For C and c in the above corollary, take Mi > Mq such that Ce 1 < 4. Once wc fix M > Mi, 
then for I is large, + Ce~ cM3/2 < 1, and hence by Q above, \T,j>i logiV^t)- 1 -21ogP(x)| < 1. Using 

\e x - 1| < (e-l)|x| for |x| < 1, 



-* 2 D t (t) - F{x) 2 \ = F(x) 2 |efc.>' log ^ (trl - 21ogFM ) 



< (e- l)F(x) 1 



^logiV J (t)- 1 -21ogF(x) 



(7.10) 



But from ( p7TT| ) and F(x) is bounded for x e K. Hence using ( |7.l| ), we obtain the result for (i) with 

new constants C, c and C(M). For (ii), we note that F(x)E(x) and F(x)E(x)~ 1 are bounded for x € M from 
(PH)-(Pl). □ 



From Proposition p.2| , Corollary p.3\ and Theorems |4.1| , [4. 2| , this Corollary immediately yields the following 
asymptotics for Poisson generating functions. 



Proposition 7.3. Let 2t = I - x(//2) 1 / 3 where x is fixed. As I — > oo, /or eac/i /ixed a, /3, 

ff(f;a) -> F 4 (x), 0<a<l, 
iffal) - Fx(x), 

Pf(f;a) -»• 0, a>l, 

P ; s (t;/3) Fx(x), /3>0. 

Le£ 4t = / - x(2Z) 1 / 3 w/iere x is fixed. As I — * oo, /or eac/i /ixed a, /3, 

if(t;a,/3) -» F 2 (x), < a < 1, /? > 0, 

ff(t;l,/3) - Fi(x) 2 , /3>0, 

P^(t;a,j3) 0, a>l, /3 > 0. 



(7.11) 
(7.12) 
(7.13) 
(7.14) 

(7.15) 
(7.16) 
(7.17) 
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Similarly, using Proposition 5.4 and Corollary 5.5 and Theorem 4.1, we have: 
Theorem 7.4. Let 2t = I - x(l/2) 1 / 3 where x is fixed. As I — > oo, we have for any fixed w € 

2^w 



if (f ; a) -F" (*;«;), 



11/3 



Let At = I — x(2t)^' 3 where x is fixed. As I — > oo, we have for each fixed f3 and aiel, 



2 5/3 

a = 1 - — tt7t — i /3 > 0. 



Jl/3 



(7.18) 



(7.19) 



Recall the relation between Qf (A) and Pf (t) in (|4.7|)-(|4.9|). We now use the de-Poissonization Lemma rcL 



to obtain the asymptotic results of Theorem |3.1| and |3.2[ In order to apply the de-Poissonization Lemma, we 
need the following monotonicity results. 

Lemma 7.5 (Monotonicity). For any I, Pr(L^ m < I), Pr(L® m < I) and Pr(L® m m _ < I) are monotone 
decreasing in k, m, m + and m_. 

Proof. We first consider Pv(L^ m < I). Let fkm '■= P r (£f m < ' |<Sfc,m| be the number of elements in Sk, m with 
no increasing subsequence greater than I. Consider the map h : Su.m-i x {1, 2, • • • , 2k + m} — > 5fc !?n defined as 
follows: for (n,j) 6 Sk,m-i x {1j2,-" , 2fc + m}, set h(n,j)(x) — w(x) for 1 < x < j — 1, h(n 7 j)(j) — j, and 
h(ir,j)(x) — ir(x — 1) for j < x < 2k + m. Then it is easy to see that /i _1 (cx) consists of m elements, hence 
(2k + m)\Sk,m-i\ — m \Sk,m\- Moreover if ir e Sk, m -i has an increasing subsequence of length greater than I, 
then h(ir,j) has an increasing subsequence of length greater than Thus (2k + m)fk<m-i) > mfkm- But since 
|5 fc , m | = f±^, we obtain Pr^,, < Z) > Pr(i^ m < I). 



A similar argument works for the other cases. Note that \S\ 



771-j- I 



(2fc+m + ,m„)! 
fe!m_j- !m_ ! 



□ 



Thus Lemma (6.3) can be applied to obtain the asymptotics results in Theorems 3.1 and 3.2. The proofs 
are similar to that in Section 9 of ||. 

Now we consider convergence of moments. For this, we first obtain the following estimates which follow 



from Proposition 5.1 (i),(ii),(iv),(v) above. The proof is very similar to that of Lemma 7.1 (i),(ii),(iv),(v) of 
Compare the results with (2.11)-(2.14) noting Corollary |7.1[ 

Corollary 7.6. Set 



t = l--l 1 /\ 
2 



(7.20) 

There exits M 2 such that for a fixed M > M 2; there are positive constants C = C(M) and c — c(M) such that 
the following results hold. 

(i). For x > M, 



c\.r\ 



3/2 



3>l 3>l 

l-l[(l-7r 2j+2 (0-t)), l-rj(l+7T 2i+1 (0;t)) < Ce~^'\ 

3>l 3>l 

l-n( 1 + ^+ 2 (0;t)), l-JJ(l-^+i(0;*)) < Ce-^'\ 



(7.21) 
(7.22) 
(7.23) 



3>l 



3>l 
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(ii). For x < —M, 



3>l 3>l 

11(1-^+2(0;*)), H(i + n 2]+1 (0;t)) < Ce- C N 

3>l 3>l 

n( 1 + ^+ 2 (0;t)), n( 1 -^ J +i(0;t)) < Ce+^l-l 

3>l 3>l 



3/2 



(7.24) 
(7.25) 
(7.26) 



Remark. From the definitions of P®(t) and the equalities of Theorem 4.1, we know that all the infinite products 
above are between and 1. 



Now as in Section 9 of p[, using Lemmap.4 and Theorems 3.1, 3.5, this implies Theorem 3.3 



Theorem 3.5 follows from Theorem 4.3 



8 Proofs of Theorems 13.41 and 13TB 



In this section, we prove Theorem 3.4 by summing up the asymptotic results of Theorems BTl p.2| and 3.3 



Theorem 3J3 can be proved in a similar way from Theorem O 
Proof of ( [3.20 ). Note that we have a disjoint union 



2k-\-m—n 



(8.1) 



Set pf ml = Pr(i^ m < I), the probability that the length of the longest decreasing subsequence of it 6 Sk, m is 
less than or equal to I. As the first row and the first column of it in S n have the same statistics, we have 



Note that 



Pr(£„<0 = 7^ T £ P f rnl \S k 



2k + m\ (2k)\ 
2k J 2 k k\ ' 



As n — ► oo (see pp. 66-67 of |p3|] ), we have 

\~S n \= \S k .^ = ^n^e-^+^-^(l + ^n-^ + 0(n-^)), 



and the main contribution to the sum comes from \fn — n e+1 / 4 < m < \fri + n e+1 / 4 . 



(8.2) 



(8.3) 



(8.4) 



Fix < a < 1 < b. We split the sum in fl8.2j ) into two pieces : 

1 



Pr(L n < 



\S„ 



(*) (**) 



(8.5) 
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where (*) is the region a\fn <ni< b^pn and (**) is the rest. 

For 2k + m = n, the quantity |5fe, m | = (2k) 1S uinm0 dal for < k < n, and the maximum is achieved 
when k ~ (n — \fn)j1 as n — > 00. Hence 



I <n-max(|5 /fe)toV s|| ) l^.^l). 

(**) 



Using Stirling's formula for (3.3), for any fixed c, when 2fc 4- [c\/n] = n, 



i^[ CV ^i~«" /2 ^" /2+ ^ (Mlogc) 



,1/4 



Hence using ([O]), we have 
1 



(8.6) 



(8.7) 



(8.8) 



(**) 



But /(x) = a; — 1 — xlogcc is increasing in < x < 1, is decreasing in x > 1, and /(l) = 0. Therefore there are 
positive constants C and c such that for large n, 



|S 



(8.9) 



On the other hand, Lemma |6.3j says that (recall (4.8)) for any fixed real number d > 0, there is a constant 
C such that for a^/n < m < b^/n 1 



(8.10) 



< P f ml <if ((2^)1/2.^(2^))-!^ +Cn - d/ ^ 
for sufficiently large n. Since Pj S (t;/3) < P^fcp), Theorem p] for P|" +1 (*;/3) yields 

^^-i)/^^) 172 ) - Cn ~ d/2 ^ P™l < e-^^f / + 2] ((2^ ) ) 1 / 2 ) + Cn-^. 
Let Z = [2y/n + iti 1 ' 6 ]. For 001 < m < b^/n, hence for "-^V" < k < "-°V" ; 

(1/2 - (2^ £i) ) 1 / 2 )2(//2)- 1 / 3 , (Z/2 - (2^) 1 / 2 )2(Z/2)- 1 / 3 = z + O^ 1 / 6 V^) 
Also note that from the asymptotics (fD|), and ( p.ll[ )-( pT^ ), 

= --(«(&) + 7x(ar))F(x)J5(x) 
is bounded for a; £ K. Hence using Q in Corollary [7^, ( pL2] ) and ( p^ ), we obtain 
|e-^D++((^)V>)-(F£0(*)| 



(8.11) 



(8.12) 



(8.13) 



< 



e-"» 2?++j((2^)V2) - (F^)((//2- (2^)i/ 2 )2(Z/2)-i/ 3 )| 

+ |(^((i/2 - K d) ) 1/2 )2(Z/2r 1/3 ) - 
< C(M)n- 1/6 + Ce~ cM3/2 + Cn~ 1/6 ylogn. 



(8.14) 
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Therefore we have 



T,Pn m i\Sn, m \ < (f(x)E(x) + C(M)n-^ 6 + Ce~ cMV2 + Cn" 1 / 6 ^^) E ( 8 ' 15 ) 

(*) V ' (*) 



Similarly, 



Y,Pnmi\Sn, m \ > f F(x)E(x) - C(M)n~ 1 ^ 6 - Ce- cM '-' - Crr 1 / 6 ^^] E \S n<n \ 
(*) V ' (*) 



But from (8.9), 



(*) (**) 
Thus using Q, dj|, ( |8T5| ), ( gig ) and ( gig ), we obtain ( |3~20| ) 



(8.16) 



.17) 



□ 



Proof of ( [3.22 ). As in Section 9 in 0], integrating by parts, 



E(fc) f 



/"Co 

pa; p - 1 ^„(x)dx+ / px p - 1 (l- J F 1 „(x))dx, 

-oo JO 



where F n (x) := Pr(x„ < ^) = Pr(i„ < 2y / n + iii 1 ' 6 ). From Theorem 4.1 and Corollary 7.6, we have 



1- e- tA ' 2 D+ + {t) < Ce~ c ^ 3/ \ x>M, 



/ 2 n++ 



(8.18) 



(8.19) 
(8.20) 



Dt + {t) < Ce- C ^ , x < -M, 

for a fixed M > M 2 where t=l - (x/2)l^ 3 . Noting that P 2 ^ +1 (i; (3) = ' /2 D+ + (t) for all > 0, from (fj 
Lemma |1| and ( ^19| ), ( |8~20| ), we obtain 



-rj.rl 3 / 2 



1-F n (x) < Ce~ cM 



x > M, 



F„(i) < Ce"^ 3 , 2<-M. 



Now using convergence in distribution, the dominated convergence theorem gives (3.22). 
Remark. We could also proceed using 

Pr(L„</) = ^ Pkml\Sk,m\. 



.21) 
.22) 

□ 

.23) 



2/c+m— n 

The main contribution to the sum from |Sfc ]m | comes from the region \m — ^fn\ < n 1 / 4+e . On the other 



hand, from Theorem 3.2, when m = s/n — 2wn 1 / 3 , the quantity p^ m [ converges to F(x;w). Since the region 
m = \/n + cn 1 / 4+€ is much narrower than the region m = y/n + cn 1 / 3 , the main contribution to the sum comes 
from when w = 0, implying that 



1 " 

Pr(L„ < I) ~ Ffa 0)\Sn,m\ = Fi(x) 



In the following proof for signed involutions, we make this argument rigorous. 



.24) 
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Proof of ( [3.21 ). We have a disjoint union 



Hence again 



One can check that 



Hence, we have 



I n I 2fc+m + +m_ — n 



| /c,m_|_ ,m_ | 



(2£ + m + + m_)! 
fc!m + !m_! 



where 



| I _ \ | 

l°n | — 2-~/ I * 

2fc-f +m_ — n 



f(m+,k) 



£ £ /(m + ,fc) 

0<fe<[§ ] 0<m + <n-2fc 



(8.25) 



.26) 



.27) 



.28) 



(8.29) 



m_)_!(n — 77i + — 2fc)!fc! 

For fixed < k < f(m+,k) is unimodal in m + and achieves its maximum when m + ~ n/2 — k. And 
/(n/2 — fc, fc) is unimodal in fc and the maximum is attained when k ~ n/2 — yJn/2. Hence /(m+, fc) has its 
maximum when (m + ,k) ~ (a/^j f — \/f )• Consider the disc -D of radius n 1 / 4+c centered at (y^, § — y/^)- 



We will show that the main contribution to the sum in ( 8.28 ) comes from D. Set 



By Stirling's formula, 



/( m fc) = ^_(2n)"/ 2 e-"/ 2+ ^ e (i + („-i/4+3e/2 )) _ 

V en7r 



;( /2 e -n/2+ v / 2H e -5v / 2n 



Hence from the unimodality discussed above, 

En 2 
f(m+,k)<n 2 max f( m+ ,k)<-^(2 
(m + ,k)£dD y/enw 



(m + ,k)<£D 



and by summing up using ( 8.31 ), 

f(m+,k) = -^(2n)"/ 2 e-"/ 2 +^ ;r (l + 0(n- 1 / 4 + 3£ / 2 )). 

(m + ,k)eD ^ 2e 

Hence we have 

\Sf\ = ^(2n)"/ 2 e-"/ 2 +^ 2;r (l + 0(n- 1 / 4+3e / 2 )), 
V2e 



(8.30) 



.31) 



.32) 



.33) 



(8.34) 
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and the main contribution to the sum in (S.28) comes from D. 



As in Theorem 3.4, we write 



Pr(L® < J) 



D" 



From (|8.32|) and ( |8.34| ), 



IcKll / fPnm+m—ll J n,m + ,ro_ 
| n | £)c 



< e 



On the other hand, by the remark to Lemma 6.3 and Theorem 4.1, (recall (4.9)) 



Pkr, 



,{d) 



(8.35) 



(8.36) 



(8.37) 



for large n. We have a similar inequality of the other direction with u, I and +Cn d l 2 replaced by fi, I — 1 and 
-Cn" d / 2 . 

Let I = [2y/2n + x2 2 / 3 (2n) 1 / 6 }. In the region D, 



and 



{1/2 - 4(^ d) ) 1 / 2 )(Z/4)- 1 / 3 = a: + O^ 1 / 6 y 7 ^), 



(1 - ^j/(4 d) ) 1/2 )2- 4/3 (V2) 1/3 = O^ 1 ^). 



(8.38) 



(8.39) 



Hence as in ( 8.14 ), using Corollary 1J2 (7.5), Proposition 5^ and Corollary 5J3 



'[1/2] 



id) 



(8.40) 



< Cn- 1 / 6 Vlo^+ Cn- 1 / 12 ^ + C(A/)n~ 1//6 + Ce" cM3/2 



for a constant C(M) which may depend on M and constants C and c which is independent of M. Thus for 
large n, 

'L*-2y/2n 



Pr 



< x < F M (x;0) + e(n,M) 



(8.41) 



2 2/3( 2n )l/6 

with some error e(n, M) such that lirnn/— >oo linin^oo e(n, M) — 0. Similarly we have an inequality for the other 
direction. Recalling F®(x; 0) = F^x) 2 from we obtain fl3.2l|). □ 



Proof of ( [3.23 ). Integrating by parts, 

/>00 

H(Xn) P )= x p dF n (x) = - I px p - 1 F n (x)dx+ I p X v- l {l-F n {x))dx, (8.42) 



where F n (x) := Pr(x^ < x) = Pr(L® < 2a/2ti + x2 2 / 3 (2n) 1 / 6 ). Note that when x < -(An) 1 / 3 , F n (x) = 0, and 
when x > 2 1 / 6 n 5 / 6 - (4m) 1 / 3 , F n (x) = 1. 
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Let M > M fixed. Consider the case when - (An) 1 / 3 <x< —M. From (|8.35|) and fl8.36|) 



n I Q K 



7714., m„ 



Ce" 



(8.43) 



where I = [2\/2n + x2 2 / 3 (2n) 1 / 6 ]. We apply Lemma Q Corollary [t7^ and Corollary |57t| fl5.36| ). Note that we 
are in the region a — ► 1 faster than fc -1 / 3 , hence w is bounded, say —1 < w < 1. So we can apply Corollary 5.7 



( 5.36 ). Then we obtain 



F n (a;) < Ce~ c W 3 + Ce 



Since — (4n) 1 / 3 < x < —M, we have 



-M 6 



thus, 



^(iJ^Ce-'W' + Ce-^W 8 ". 
On the other hand, when M < x < 2 1 / 6 n 5 / 6 — ^n) 1 / 3 , similarly we have 



1 - F n (x) < T^^O—Pkm+m-OlSf^ 
o, 



+ Ce 



n I D 



Using Lemma 5A, Corollary 'Lt and Proposition 5^ ( 5.32| ), we obtain 

l-F n {x)<Ce- c ^ m +Ce- Wn * 

Since M < x < 2 1 / 6 n 5 / 6 - (4n) 1 / 3 , 



thus 



< e 2 2 



1 - F n (x) < Ce-°l x l 3/2 + Ce~^ 



i° \x\ 12 '^ 



10 |_|12e/B 



Therefore, using dominated convergence theorem, we obtain (|3.23| ) 



(8.44) 



(8.45) 



(8.46) 



.47) 



(8.48) 



.50) 
□ 



9 Asymptotics for a > 1 



As we remarked after Theorem 3.3, when a > 1, we must use a different scaling to obtain useful results. 

Let L (i;a) and L®(t;a,(3) be random variables with the distribution functions given by Pr(L (i;a) < 
I) = P (t; a) and Pr(L H (t; a, 0)<l) = Pf (t; a, /3), respectively: the Poissonized version of L and L . Under 
appropriate scalings, we obtain the Gaussian distribution in the limit. 
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Theorem 9.1. For a > 1 and (3 > fixed, 



f L a (t;a) - (a + a~ l )t \ 1 f x i„a , , , 

lim Pr v ' / v = — i- < x = -= / e~2 y dy, 9.1) 



X 



,. f L ft; a, f3) — 2(a + a )t \ 1 , 

lim Pr V ; ^ V ^ L < x = __ / e ~2V dy. 9.2) 



Proof. Let 2 = (a + a x )i + ^/{a — a~ x )t for L . For large i, 2t/2 < c < 1 for some c > 0. From Proposition 
5.1 (i), using Theorem 4.1, it is easy to see that e~ l l 2 D^(t), e~ l ± *D ; ±=F (t) — > 1 exponentially as 2 — > oo. Now 



Theorem gj, (|5.29|) and Proposition |J imply (|9T|) . For L H , let 2 = 2 (a + a -1 )* + y / 2(a - oF 1 )!. Similarly, 
e~* Di(t) — ► 1 exponentially as 2 — * oo, and we obtain (|9.2[). □ 



Unfortunately, we can no longer apply the de-Poissonization technique; the difficulty is that (a + a _1 )< 
depends too strongly on small perturbations in a. Indeed, as we shall see, the asymptotics of the non-Poisson 
processes are different. 

Consider the case of involutions with [at] fixed points and [t 2 /2] 2-cycles; the case of signed involutions is 
analogous. By symmetry, this is the same as the largest increasing subset distribution for the point selection 
process in the triangle < y < x < 1 with [t 2 /2] generic points and [at] diagonal points. As was observed in 
Remark 2 to Corollary 7.6 of 0, it is equivalent to consider weakly increasing subsets where the extra points 
are added to the line y = instead of to the diagonal. 

As in fllTi"!) , let 

v0 _ L \tV2],wt] - ( Q + W rq „x 

X[t2/2] < [Qt] - vd/«-v« 3 )i ■ ( 3 

Theorem 9.2. ^4s £ — > oo, the variable X^/ 2 ] [ a t] converges in distribution and moments to N(0, 1). 

Proof. Let S(t) be the set of points at time t, and let / be a largest increasing subset of S(t). Then there will 
exist some number < s + < 1 (not unique) such that 

(S(t) n {y = 0, < x < s+}) C I, (9.4) 

and such that every other point of I has x > s + and y > 0. For any < s < 1, we thus have 

fi(s) + f 2 (s)<\I\, (9.5) 

where /i(s) is the number of points of S(t) with y = and < x < s, and where /2(a) is the largest increasing 
subset of S(t) lying entirely in the (part-open) trapezoid with x > s, y > 0. 
Since /i(s) is binomial with parameters [at] and s, 

Lemma 9.3. Let M > be sufficiently large and fixed. For all < s < 1, there exist positive constants C, c 
independent of s such that for w > M, 

Pr(/i(s) > asi + wt 1/2 ) < Ce~ clwl \ (9.6) 

while for w < — M, 

Pr(/i(s) < as< + urf 1/2 ) < Ce-^"" 12 . (9.7) 
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For /2, we have: 

Lemma 9.4. Let M > be sufficiently large and fixed. For < s < 1, there exist positive constants C and c 
independent of s such that for all w > M , 

Pr(/ 2 (s) > 2V(1 — ~s)t + wt 1 ' 3 ) < Ce~ c ^ 3/ \ (9.8) 

and for all w < —M, 



Pr(/ 2 (s) < 2^(1 -s)t + wt 1 ' 3 ) < Ce- c|m| . (9.9) 

Proof. We first show the corresponding large-deviation result for the Poissonization. Define /^(Sji) to be 
the length of the longest increasing subsequence when the number of points in the trapezoid is Poisson with 
parameter i 2 (l—s 2 )/2. Then f^is , t) is bounded between the corresponding processes for the rectangle s < x < 1, 
< y < 1 and for the triangle < (x — s)/(l — s) < y < 1. In particular, if f' 2 (s,t) deviates significantly 



from \/l — st, so must the appropriate bounding process; the result follows immediately from the corresponding 
results for rectangles and triangles. 



The corresponding large-deviation result when the number of points is fixed then follows from Lemma x2 
In our case, the number of points in the trapezoid is binomial with parameters t 2 jl and (1 — s 2 ); the lemma 



follows via essentially the same argument used to prove Lemma S.2. □ 

As we will see, the value s = 1 — or 2 deserves special attention: 
Lemma 9.5. Let M > be sufficiently large and fixed. There exist positive constants C, c such that for w > M, 

Pr(/x(l - a' 2 ) + f 2 (l - a' 2 ) > (a + l/a)t + wt 1 ' 2 ) < Ce~ cmiD tt w \ 2 ^>l 3/2 * 1/4 ), (9.10) 

and for w < — M , 

Pr(/i(l - a~ 2 ) + / 2 (1 - a- 2 ) < (a + l/a)t + wt 1 ' 2 ) < Ce^l 2 . (9.11) 

Moreover, if we define 

(t) = Ml-^) + h(l-^)-(a + l/a)t 
v/(l/«-l/a 3 )i 

then Xo(t) converges to a standard normal distribution, both in distribution and in moments. 

Proof. That xo{t) converges as stated follows from the fact that if we write Xo(£) = Xi(^) + with 

() _h(l-a- 2 )-(a-l/a)t 
^(l/a-l/a 3 )t 

X 2 (t)= hil - a ' 2) - {2/a) \ (9.14) 
' y/(l/a- l/a 3 )t K ' 

then xi{t) converges in distribution and moments to a standard normal distribution, and x.2{t) converges in 
distribution and moments to 0. 



3G 



For the large deviation bounds, we note that if x + y > z + w, then cither x > z or y > w. Thus for any 
< b < 1, we have 

Pr(/i(l - a- 2 ) + / 2 (1 - a- 2 ) > (a + l/a)t + wt 1/2 ) (9.15) 

< Pr(/i(l - a- 2 ) > {a- l/a)t + bwt 1/2 ) + Pr(/ 2 (1 - a" 2 ) > (2/a)t + (1 - b)wt 1/2 ) (9.16) 
< Ce -c|6H 2 + C - e -c|(i-6) w |W4 ; (91?) 

the result follows by balancing the two terms. In the other case, the Ce _c ' t0 ' 2 term always dominates. □ 

Lemma 9.6. For any sufficiently small e > 0, there exist positive constants C , c such that 

Pr(s+ - (1 - l/a 2 ) > ^/ 3 -V3) < Ce~ ct ' (9.18) 

and 

Pr((l - l/a 2 ) - s+ > f' 2 - 1 ' 2 ) < Ce~ ce (9.19) 

for all sufficiently large t. 

Proof. Define a sequence Si by taking 

s t = l-{l-2/{i + 2)) 2 /a 2 (9.20) 
for all i > 0. Similarly define a sequence s- by 

s'i = max(ti,0), (9.21) 

with 

U = 1- (l + 2e 2 ~ 1 ~') 2 /a 2 (9.22) 

for i < and 

k = l- (l + 4/(t + l)) 2 /a 2 (9.23) 

for i > 0. Note that Sj is strictly decreasing and ti is strictly increasing. 
Lemma 9.7. For all i > 0, 



as, + 2^J\ - s i+1 <a + l/a. (9.24) 

For all i, 

as' l+1 + 2^J\ - s\ <a + l/a. (9.25) 
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Proof. In the first case, we have 



a + l/a- (asi + 2y/l - s l+1 ) 



{i + 2) 2 {i + 3)a 

In the second case, it suffices to verify the formula with s' replaced by t. For i < — 1, 

a + l/a - (at l+1 + 2^/l~~t~) = 4e 2 ~ V4 /a. 

For i = — 1, 



a + l/a - (at l+1 + 2\/T — 1~) = (24 - Ae)/a 

Finally, for i > 0, 

8i 



a + l/a- {aU+i + 2Vl~U) = 



{i + l){i + 2) 2 a 



Let ii — t 1 / 6 e / 6 , %2 — t 1 / 4 e / 4 . Then there exist constants C, c such that for < i < i\, 
Pr(/i(«i) + / 2 («i+i) > /i(l - O + / 2 (1 - a- 2 )) < Ce~ ct \ 
Since /i(si) + /2(si+i) is an upper bound on /i(s) + /2(s) with Sj+i < s < Sj, it follows that 

Pr(s+ e [s i+ i, Si ])<Ce- c * e . 

Similarly, for z < i 2 , 

?*(fM+i) + Ms-) > /i(l - a" 2 ) + / 2 (1 - a- 2 )) < Ce~ ct \ 

and thus 

Pr(s + el S >, S > +1 ])<Ce- ct '. 
Since there are only i\ + i 2 + log log a such events to consider, the result follows. 
In particular, with probability 1 — Ce~ ct ' , we have 

/i(l - a- 2 - f-/ 2 - 1 ' 2 ) + f 2 (l - a- 2 + t«/ 3 -V3) < L(t) < A(l - ^ 2 + W 3 - 1 / 3 ) + / 2 (1 - a" 
But then using the fact that 

/l (l-a- 2 )-/ 1 (l-a- 2 -^/ 2 -V2 ) 

and 

/i(l-«- 2 -t £/2 - 1/2 )-/i(l-«- 2 ) 
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are Poisson, and using the large deviation behavior of /2(s), we find that 

Pr(L(t) - - a" 2 - i £ / 2 -V2) + / 2 (! _ a -a + ^/a-i/ajj > t Va- £ ') < c^-of j (9.37) 

and 

Pr((/i(l - a" 2 - ^-VS) + ; 2 (i _ a -J> + ^/a-i/ajj _ L(i) > fi/a-e'j < Ce -ct« (9 38) 

So x(i) — Xo(^) converges to in a fairly strong sense; in particular, they must have the same limiting distribution 
and limiting moments. □ 

Remark. The above proof could equally well be applied to the Poisson process; the beta distribution would 
then be replaced by a Poisson distribution. 

For the signed involution case, with [2at] fixed points, [2(3t] negated points, and [2t 2 ] 2-cycles, again we let 

X[«"].M.PW " V 2(l/a~l/a*)t ' ^ 

Then the analogous argument proves that x?{t) also converges in moments and distribution to a standard 
normal. 

10 Steepest descent type analysis for Riemann-Hilbert problems 

In this section, we prove the asymptotics of orthogonal polynomials results stated in Section [| by applying the 
steepest descent type method to the Riemann-Hilbert problem (RHP) ( |5.4| ). The steepest descent method for 
RHP's, the Deift-Zhou method, was introduced by Deift and Zhou in Jl8| and is developed further in JlOf] , |lg| , 
and finally placed in a systematic form by Deift, Venakides and Zhou in [jl7|. The steepest descent analysis of 
the RHP (5.4) was first conducted in Q. The analysis of || has many similarities with [jl3|, Q [l5| where the 



asymptotics of orthogonal polynomials on the real line with respect to a general weight is obtained, leading to 
a proof of universality conjectures in random matrix theory. 

As mentioned in the Introduction and Section |^, in this section we extend the analysis of || and obtain new 
estimates on the orthogonal polynomials %k{z;t). The extension is done roughly in two categories. In ||, the 
quantity of interest was l2i(0; k; t), and so the z-dependence of the error bound of Y(z; k; t) was not considered 
carefully. But in the present paper we need the asymptotics of Y(z) for general z E C, and also for the case 
when z — > — lasfc,i— >oo. Hence the first category of our extension is to investigate how the error estimate 



depends on z. This task sometimes requires improved estimates of the solution Y(z). See for example, (10.42) 
below where an improved L 1 norm bound of the jump matrix is needed. On the other hand, as we will see, the 
asymptotic solution Y(z) is expressed in terms of the so-called the g-function. Thus we need detailed analysis 
of the g-function to obtain the asymptotics of the orthogonal polynomials. In the special case z = 0, we have 
g(Q) — iri (Lemma 4.2 in S). Hence in ||, the analysis of the ^-function was quite simple. But in the present 
paper, we need general values of g(z), and this in some cases requires further analysis. Hence the analysis of 
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the g-function is the second category of our extension. For example, see (10.109)-(10.121) below where we need 
a further analysis of the ^-function. 

Again the analysis in this section relies heavily on the analysis of || and we extend the method of jj| . Thus 
the analysis below overlaps in many parts with the analysis of However, for continuity of presentation and 
also for the convenience of readers, we include some calculations which overlap When the analysis overlaps 
that of ||, we only sketch the method, and instead we focus on new features to indicate how to prove the 
propositions in Section pL 



We say that a RHP is normalized at oo if the solution m satisfies the condition m — ► I as z 



Thus for 



instance, the RHP's ( 2.15 ), ( 10.1 ) are normalized at oo, while the RHP (5A) is not. 

In H), it turned out that the asymptotic analysis differs critically when (2t)/k < 1 and (2t)/k > 1, due to the 
difference of (the support of) the associated equilibrium measure (see Lemma 4.3 j})). Hence we discuss those 
two cases separately in Subsections 10.1 and 10.2, which extend Section 5 and Section 6 of j3|, respectively. 
Each section is also divided into three subcases. In each subcase the corresponding case of the propositions in 
Section || (except Proposition pTsj ) are proved. Subsection 10.3 is new, and Proposition 5.8 is proved in this 
subsection. 



10.1 When (2t)/k < 1. 



The following algebraic transformations ( lO.l )-( lO.lO ) of RHP's are taken from (5.1)-(5.3) of 
Define 



m (1) (z;fc;f) := Y(z;k;t) 



'{-l) k e tz 
. 







k ^—tz 



i {1) {z;k-t) := Y(z;k;t) 



k„tz 



z e 

. o 



\z\ <1, 



1*1 > I- 



(10.1) 



Then solves a new RHP which is equivalent to the RHP (5.4) in the sense that a solution of one RHP 
yields algebraically a solution of the other RHP, and vice versa: 



(z; fc; t) is analytic in C \ E, 

^_^k z k e t{z-z- 1 ) 

(_ 1 )fe z -fe e -t( 2 - z - 1 )^ 
m w (z;M) =1 + °{h) asz-^oo, 



m^P (z;k;t) = m^' (z; fc; t) 



on S, 



(10.2) 



where £ is the unit circle oriented counterclockwise as before. Here and in the sequel, m+(z) (resp. m_(z)) is 
understood as the limit from the left (reps, right) side of the contour as one goes along the orientation of the 
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contour. Now we define (z; k; t) in terms of m^ 1 ' (z; k; t) as follows: 

for even k, 

m (2 ) = m (1 ) \z\ > 1, 

mW = m( 1 )(?- 1 ) 1*1 <!■ 
for odd fc, 

mP) = (J_ 1 )m( 1 )(J_° 1 ) W>1, 

Then m( 2 )(-; fc; t) solves another RHP 

m?^; k; t) = m^\z; k;t)v^ 2 \z; k;t) on E, 
i (2 \z;k;t) =J + 0(i) as z^oo 

where 

(_i \k k t(z-z~ v ) 



(10.3) 



(10.4) 



The jump matrix has the following factorization, 

^(-l)^-fee-^-^ 1 ) l/ [p l y v-; +, v 7 

We note that through the changes Y — > — > m^ 2 \ we have 

Yn(z;fc;i) = -(-l) k e- tx m${z;k;t), \z\ < 1, (10.7) 

Fai («;*;*) = -e- ta m^ («;*;*), l«l<i. (10.8) 

yn(z;fc;t) = z'e-^mli'^fc;^ |*| > 1, (10.9) 

Y 21 (z;k;t) = (-a)^ -1 *" 1 ^? («;*;*), > 1. (10.10) 

As in (5.4) of ||], the absolute value of the (12)-entry of the jump matrix is e kFt > p,e '^ where 

F(z; 1 )=F(pe ie - n ):=l(p-p- 1 ) cos6 + \ogp, z = pe w . (10.11) 

The absolute value of the (21)-entry of i/ 2 ) is e - kF (P e *' Note that 

F( P ,0;j) =-F{p-\6; 1 ). (10.12) 

Figure | shows the curves F{z\i) = 0. In ftf 9 U fl s 3 w , F > 0, and in fi^ 9 U Q s 4 19 , F < 0. The region ft^ 9 
becomes smaller as 7 increases, and when 7 = 1, the curve F(z;j) = contacts the unit circle E at z = — 1 
with the angle 7r/3. 

We distinguish three subcases as in Section 5 of §. 
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Figure 3: curves of F(z; 7) = when < 7 < 1 



10.1.1 The case < 2t < ak for some < a < 1. 

It is possible to fix p a < 1 such that the circle {z : \z\ — p a } is in the region Slff 9 f° r an such t and k. Define 
m (3 \z;k;t) by (see (5.9) of @) 

' m (3) =m (2)( 6 W)-l j p Q <|z|<l, 

m (3) =m (2)( 6 (2))-l ) 1< | Z | < p- 1 , (10.13) 

m( 3 )=m( 2 ) kl<Pa, M > Pa 1 - 

Then m^ 3 ) satisfies a new jump condition m^ 5 = mf 3 '?/ 3 ) on X^ 3 ) := {z : |z| = pa,/^ 1 }, where = b+\ 
\z\ = p a and z/ 3 ) = (&_ \z\ = p" 1 . This E^ 3 * 1 is not the best choice (see the next subsubsection). But 
for a simple and direct estimate, we use this choice in this subsubsection. From the choice of p a , we have (see 
(5.13)^(5.14) of |) 

|v (3) (z;fc;i) -I\ < e~ ck for all z G £ (3) , (10.14) 

which implies that I—C w (a) is invertible and the norm of the inverse is uniformly bounded, where := i/ 3 -* — I 
and C wi3) (f) := C-(fw^>) on L 2 (£< 3 ), \dz\), C± being Cauchy operators (see (2.5)-(2.9) in H and references 
therein). From the general theory of RHP's, we have 

2m J S (3) s- z 

This implies the estimates (see (5.16) of 0|) 

\m^(0;k;t)-l\, |mg } (0; k; t)\ < Ce- ck , (10.16) 
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which using (10.13), (10.7), (10.8) and (|5.5|), (5.6), yield Proposition 5.1 (i). This is precisely the result contained 
in (5.17) of §. 



From (10.13), (10.7), (10.9) and (§J), we have 



,(3)/ 



ir k (z;t) = 


-(-1) 


TTk(z;t) = 


z k e~ tz 


ir k (z;t) = 


z k e~ tz 



\z\ < Pa, 



(10.17) 
(10.18) 
(10.19) 



Let < b < 1 be a fixed number. From Figure |[ we could have chosen p a such that p a > b. When \z\ < b and 
\z\ > we have dist(z, E^ 3 -*) > c > 0. Since z is uniformly bounded away from the contour, we can extend 
the argument leading to (5.17) in [pi where the uniform boundedness of from the contour is used. Hence using 
(|10.15p, (|10.17|) and (|10.18j) imply that 



\e tz n k (z;t)\ < Ce- ck , 
'z- k ir k (z;t)\ < Ce- Ck , 



\z\<b, 
\z\>b-\ 



(10.20) 
(10.21) 



These are ( 5.29| ), ( ^5 .30 ) in Proposition 5.6 of the special case x > 2 1 / 3 (1 — a)k 2 ^ 3 . 

On the other hand, let L > be a fixed number. Set a = 1 — 2 4 / 3 fc~ 1 / 3 ui with — L < w < L as in Proposition 
5.6 . Since p a is fixed, when k is large, dist(— a, X^ 3 )) > c > 0. Then from ( 10.15 ), 

\m$(-a;k;t)-l\, \m%\-a;k;t)\ <Ce~ ck . 



Note that 



-(s-s- 1 ) < s-1, s>0, 



--(s - s _1 ) +logs < 0, S>1. 



5 <-<i, 



(10.22) 

(10.23) 
(10.24) 
(10.25) 



Thus for 7 < 1, s > ± 



logs 



1-7, -In 



-X s 



logs 



<(l_ 7)(s _l ) + _| s _l|3. 

For large k, a> h for all — L < w < L, and hence 



(10.26) 



(—a) e v ~ e 'fc^<e 2 * < e 2 e 



Similarly, since a 1 > ^ for large fc, 

\(~ a )- k e^ a - a ^\ 



#L 3 2 4 / 3 Lfc 2 / 3 



= Ce 



cfc 2 / 3 



fcF(-a- 



^)</(i-f)^ + fHV» <c , 



* 2 / 3 



(10.27) 



(10.28) 
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Therefore from fltO.lSD and ( |10.22| ), 

\e~ ta n k (-a]t)\ = \{-a) k e- t{a ' a ' 1) m{\\-a]k;t)- (-l) k mf 2 \-a:k;t)\ < Ce ckV \ (10.29) 
|e-* Q_1 (-a)- fe 7r fe (-a;t) - l| = \m { ^ {-a;k;t) - 1 - a^e^-^m^ {-a;k;t)\ < Ce~ ck . (10.30) 

Noting x ~ A: 2 / 3 , these are (5.31), ( 5.32| ) in Proposition 5.6 of the special case x > 2 1 / 3 (1 — a)fc 2 / 3 . Thus we 
have extended the argument of (5.17) of || to the case when z — > —1. This is an example of the extension of 
the first category mentioned at the beginning of Section fffl (though it is straightforward to extend in this case) . 



10.1.2 The case ak<2t <k- M2~ 1 / S k^ s for some < a < 1 and M > M . 

In the previous section, the contour was not the best choice. We could have chosen the steepest descent 
curve for F(z;j). For the previous case, it was not necessary to use the steepest descent curve to obtain the 
desired results, but for the case at hand and in the future calculations, we need to use the steepest descent 
curve. 

For fixed 9 satisfying < 9 < ir/2 or (37r)/2 < 9 < 2ir, F(p,9;-f) is always negative for < p < 1, and as 
p I 0, it decreases to minus infinity. On the other hand, one can check that (see (5.5) in ||) when j < 1, the 
minimum of F(p,9\ r f), < p < 1 is attained, for fixed tt/2 < 9 < (3ir)/2, at 



1 - a/1 - 7 2 cos 2 9 



-7 cost 



and F(p e , 9; 7) < 0. And it is straightforward to check that for < 7 < 1, tt/2 < 9 < 3tt/2, 

2V2, 



F(p e ,9; 1 ) = y/1 - 7 2 cos 2 9 + log 



1 — \Jl — 7 2 cos 2 



-7 cos ( 



< -- 



-(1 + 7 cos6>) 3/2 . 



(10.31) 



(10.32) 



This is an extension of (5.13) in |jj where only the case 9 — it was considered. We need this improved version 
to obtain a better L 1 estimate of — I in the sequel. Also F(po,9;j) is increasing in tt/2 < 9 < tt and is 
decreasing in tt < 9 < 3ir/2. In fact, the saddle points for \{z — z^ 1 ) + log z are z = — and z = —p^ 1 - 



In this time, define £< 3 ) := U by, as in (5.6) of 



,(3) 



= {p e e w : 3tt/4 < 9 < 5tt/4} U {p 3v/i e w : < 9 < 3tt/4, 5tt/4 < < 2tt}, 
£?Jt = {Pe 1 ^ ■ 3tt/4 < 9 < 5tt/4} U {p^ /4 e W : < 9 < 3tt/4, 5tt/4 < 9 < 2vr}, 



(10.33) 



where p e is defined in ( |l0.3l| ) with 7 = (2t)/k. Orient E^ as in Figure |. Note that £ CT lies in fi^ and for 
37r/4 < 9 < 57r/4, it is the steepest descent curve. The reason why we choose a part of the circle as the contour 
for the remaining angles is to ensure the uniform boundedness of the Cauchy operators (see Section 5 of ||.) 
This does not affect the asymptotics since as we will see the main contribution to the asymptotics comes from 
neighborhood of z = — 1. 

Define the regions j = 1, • • • , 4 as in Figure 0. Define m^(z; k; t) by, as in (5.9) of ||, 



' m (3) =m (2)( 6 W)-l j 
m (3) =m (2)( 6 W)-l ) 

m' 3 ' = 7?i/ 2 ) 



in f2 
in Q 



(3) 

2 ' 

(3) 

3 ) 



(10.34) 



n (3) 

\ 1 J '4 ; 
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Figure 4: T,^ and Q( 3 ) when 7 < 1 



where b±^ are defined in ( 10. 6| ). Then solves a new RHP with the jump matrix v^(z; k;t) where 



,(3) 



,(3) 



1 -{-l) k z k e^ z - z *) 
,0 1 

1 0' 

(-l^z-ke-Kz-z- 1 ) j 



on S 



(3) 



(10.35) 



on £ 



(3) 



Set u/ 3 - 1 := w' 3 -' — /. For z G from the choice of E^ 3 ^ and (10.32), the (12)-entry of the jump matrix 

satisfies for 37r/4 < argz < 57r/4, 



and for < argz < 37r/4 or 57r/4 < argz < 27r, 



5 e) 3/2 



< e 



-2#fe(l-f) 3/2 < e -|M 3 / 2 



l^kgtC*-.- 1 )! = e fcF(p 3 „ /4 ,e;f ) < e fc J F(p 3W4 ,3 | L ; f ) < e _lV2 fe(1+ 2i cos ^ ) 3/2 ^ 

From ( 10.12 ). similar estimates hold for z - fc e -*0- z ), 2 € Thus we have 

lk (3) lL«., v . 



(3) a < Ce — k{1 s ^ . 



(10.36) 



(10.37) 



(10.38) 



Also there exists Mq such that for M > Mo, HC^s) \\l 3 Ce(.*))->l 2 Ce( 3 ')) < ci < 1, and hence ( 10.15 ) holds. This 
is precisely (5.18) of ||. For this derivation, we do not need the extension (10.32) of (5.13) in ||. But for the 
improved L 1 norm estimate of w^, which we do now, we need ( |10.32[ ). 

Note that \dz\ < C\d0\ on S< 3) . Using the estimates in (|10.36[ ) and ( |10.37| ), 



\ z k e t(z-z-i)\ 



\dz\<C e-W+xcos^^ + c 



e 2Z k dd. 



[0,27r)\[3ir/4,57r/4] 



(10.39) 
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The second integral is clearly less than Ce~^s . For the first integral, recall the inequality (x + y) a > x a + y a , 
x,y > 0, a > 1. Then using the inequality 1 + cos 9 > 7^75(6* — k) 2 for 8 E [=^, together with the condition 
ak < 2t, the first integral is less than or equal to 



-(f (l+cose))' 



ld»<e-W 1 - 



\ 3/2 



e 3-23/ 4 cty, 



(10.40) 



where last inequality is less than or equal to Ck for some constant C > 0. Therefore adjusting constants, 
we obtain 



\^-^\\ dz \<^e~^ 



3/2 



(10.41) 



(3) 

We have similar estimates on E^. Therefore, 



fc l- 



(10.42) 



This is a refinement of (5.23) in jjj. 
Now from (10.15), we have 



m (3) (z) = 7 + 



2m 



S(3) s — z 



2tti 



[(/-C^o))- 1 ^^/]^ 3 )^) 



ds, z^S (3) . 



(10.43) 



In ||, only the term (z) (at 2 = 0) was of interest. But then — 0, and the first integral in ( 10.43| ) was 
zero. As computed in (5.20) of ||, the second integral is bounded by the product of the L°° and L 1 norms of 
vj( 3 \ and then due to (10.38), the estimate ||iz/ 3 )||£i < Cfc -1 / 3 in (5.23) of was enough to control ■ 



,(3) 



But 



in the present paper, we need estimates of m 12 for TTk(z), and hence we need an estimate of the first integral 
which is same as a bound on the L 1 norm of u/ 3 ). The L 1 bound (5.23) obtained in |3| is not good enough for 
this purpose, and we need an improved estimate on the L 1 norm of w^ 3 \ Now by (10.42), the first integral is 
less than or equal to 



C 



3/2 



dist^E^fcVs 

while, as in (5.20) of ||, the second integral is less than or equal to, by using ( 10.38 ) and ( 10.42 ), 

1 



(10.44) 



27rdist(z,£( 3 )) 
1 



< 



< 



< 



< 



27rdist(2,£( 3 )) 
C 

dist(z,£( 3 ))" U 
C 



- c wm )- x c w ^i\\ L 4w^\\ L , 



(3) II 2 . 



(10.45) 



dist(z,S( 3 )) 

c 



distfoEO))*; 1 ^ 



lk (3) IMk (3) IU* 



4G 



When z = 0, dist(z, S (3) ) > a > 0, hence using (|10.34|), (|10.7|), (|10.8|) and (|5.5|), (pjj), we obtain Proposition 



yi(h). 



As in (10.17)-(10.19), from (10.13), (10.7), (10.9) and (5.6), we have 



nk(z;t) = 


-(-l) k e- tz m^(z;k;t), 




nk(z;t) = 


z k e-~ tz ' 1 m ( ^{z-k;t), 


zed£\ 


Kk{z;t) = 


z k e- tz ~ 1 m < ^{z;k;t) - {-l) k e- u m${z;k;t), 


z G fl^ U 



Define x by 



2t _ x 

J ~ 1 ~ 2V3jfc2/3 



(10.46) 
(10.47) 
(10.48) 



(10.49) 



as in Proposition |5.6| . Let < b < 1 be a fixed number. Given 6, from the beginning, we could have chosen 
< a < 1 such that 



P6t 



1 - VI - a 2 cos 2 b 
—a cos 9b 



(10.50) 



is strictly greater than b for some n/2 < 9b < it. Note that in ( 10.33 ) defining S^ 3 \ the choice of 37r/4 and 57r/4 
was arbitrary. Instead of 37r/4 and 57r/4, this time we use 9b and 2n — 9b, and carry this forward through the 
later calculations. Thus we obtain the same estimates of ( |10.44| ) and ( 10.45] ) with different constants C. Now 
\z\ < b lies in f^ 3) and \z\ > lies in f^ 3) . Since the distance dist(z, E< 3 >) > c 2 > 0, using (|10.43|) , ( |10.44|) , 
( |10.45[ ), ( |10.46| ) and ( |10.47| ), we obtain fl5.29| ) and (|3cj) in Proposition |J for M < £ < (1 — a)2 1 / 3 fc 2 / 3 . Since 



in ( 10.20 ) and ( 10.21 ) in Subsubsection 10.1.1 , the choice of < a < 1 was arbitrary, we obtain ( 5.2S ) and 
(5.30) in Proposition 5.6 for all x > M. 

On the other hand, let < L < 2~ 3/2 v / M be a fixed number. Set a = 1 - 2 4/3 fc~ 1/3 w with -L < w < L as 



in Proposition 5.6. From the inequality 

1 - 



< 1 — \/l — 7 for all < 7 < 1, we have 



P-K = 



2± 

k 



< 1 



2t 



M 



2 l/6 fc l/3 * 



But 



2^w 2 4 / 3 L 
a = 1 -JT- > 1 - 



jfel/3 " fcl/3 ' 

Hence dist(-a, S^ 3 )) > CAT 1 / 3 . Thus from (|l0T|), ( ^044] ) and ( p0~45| ), we have 

|m( 3 )(-a;fc;i)-/| < Ce"i* 3/2 , 



which together with (10.48|) (note that — a e Sl^ 3 '* U ^3 J,I )> implies that 



|e-* a 7r fc (-a;i)| < | (- a )* e - t(Q - Q ^ | (l + Ce-^' 2 ) + Ce~ c ^ 
(-a)- k ir k (-a;t) - 1\ < Ce^ 3 ' 2 + U-a^e^^lCe-^ 2 . 



(10.51) 



(10.52) 



(10.53) 



(10.54) 
(10.55) 
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For large k, a > | for all — L < w < L, hence using (10.26]), we obtain as in (10.27) and (10.25), 



(10.56) 



and 



|(-a)- fc e t ( a -°'" 1 )| = e*^" 01 " 1 '^) < Ce c K 



(10.57) 



Thus from (10.54) and (10.55), we obtain (5.31), (5-32) in Proposition 5.6 



10.1.3 The case k - M2~ 1 / 3 k 1 / 3 < 2t < k for some M > 0. 



In this case, as k — > oo the point z = — p n on the deformed contour E^ 3 ^ defined in (10.33) approaches z = — 1 
rapidly, and so we need to pay special attention to the neighborhood of z = — 1. More precisely, we need to 
introduce the so-called parametrix for the RHP around z = — 1, which is an approximate local solution. 






Figure 5: E p// > 2 and nf 11 ' 2 

Recall the RHP ( 2.1 5| ) for the Painleve II equation. Let E p// - 2 = E P/7 ' 2 U E P// ' 2 be a contour of the general 
shape indicated in Figure |[ Asymptotically for large z, the curves are straight lines of angle less than ir/3. See 
the paragraph after (2.18) in || for more precise discussions on the curve. We will define the exact shape of 
£ p/p2 below. Define m PII ' 2 (z;x) by 



m ' (z, x) — m(z; x) 



m ' (z, x) — m(z; x) 



m (z, x) — m(z; x) 




9. 



PII,2 



PJ/,2 







• r>PII,2 r>PII.2 
BlSlj , S2 4 , 



(10.58) 



where m(z; x) is the solution of the RHP for the PII equation given in |H|). Then m p// > 2 solves a new RHP 
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(see (2.19) of {§) 



,PII.2 



PI 1,1 PII,2 

m, = to 



to. 



PII,2 




in S 



PII.2 



(10.59) 



PII,2 



in £ 



PII,2 



i PII *=I + 0(± 



PI I 2 m RI 1 ' 2 (x) 

Also m 1 ' (x) defined by m PII - 2 (z; x) = I + m * 



0(z 2 ) satisfies 



PII.2, \ / \ 

to x (re) = toi(cc), 



(10.60) 



where mi(cc) is defined in a similar manner in (2.16). 
Set x by 



2t 



1 



2 l/3 fc 2/3- 



(10.61) 



We define S< 3 ) and m( 3 ) as in 010.3$ and (|l0.34j) . Now we proceed as in (5.25)-(5.35) of f|. Let O be the ball 
of radius e around z = — 1, where e > is a small fixed number. Define the map (see the equation displayed 
between (5.26) and (5.27) of [|) 



X(z) :=-z2- 4 / 3 fc 1 / 3 i(z-z- 1 ) 



(10.62) 



in O. Define E p// < 2 by £ P7/ ' 2 n A(0) := AfE^'nO) and extend it smoothly outside A(O) as indicated in (5.29) 
of ||. And define m PI1 ' 2 as above using this contour. Now we define the parametrix by (see the equation 
displayed between (5.30) and (5.31) of H) 



to p (z; k; t) = m PII ' 2 (X(z),x) 
m p (z; k;t) = I 

It is proved in (5.25)-(5.34) in j|] that if we take e small enough but fix it, then the ratio 



inO\£( 3 \ 
in O c \^ 3 l 



solves a new RHP 



R(z;k;t) := to^to" 1 



R(z; k; t) is analytic in C \ 

R+(z; k;t) = i?_ (z; k; t)vn(z; k; i), on S/j, 

R(z;k;t) = l + 0(i), as z -> oo, 



(10.63) 



(10.64) 



(10.65) 
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where Hr := dO U S^ 3 \ and the jump matrix satisfies (see (5.34) of 



\vr - < P73 



- ck 



\vr-I\\l~ < Ce 

I 7" i Tfh-t (X) || ^ C 

A( Z ) IU°° ^ 



on On £( 3 \ 

on o c n yP\ 

on <90, as k — > oo, 



(10.66) 



with some positive constants C and c which may depend on M. Set u># := vr — I. Using ( 10.15 ) which holds 
generally, we have (see (5.35) and the preceding calculations in Q) 

(10.67) 



R(z;k;t) = I+^- f 



1 

2ttI 



v R (s)-I 1 f [(I — C WR )~ C wr I](s)wr(s) 



-ds , 

s — z 2ni 



ds. 



Now the absolute value of the second integral is less than or equal to (recall that |A(z)| = 0(k 1 / 3 ) for z £ dO) 

C 



dist(z, Yjr) 
C 

< 



< 



dist(z, Hr) 
eC 



(I - C WR ) 1 \\l 2 (T, r )^L^{S r )\\C Wr I\\l 2 {T, r )\\wr\\l' 2 (T, r ) 

w r\\ 2 l*{-z r ) 



(10.68) 



dist(z,£ K )fc 2 / 3 ' 
and similarly, the first integral satisfies 



-/ 

2tt« Jy, r 



1 \ T PH,2/ \ 

vr{s) - I_ ds m 1 ' (x) 



27T« Jdo H s )( s - z ) 



-ds 



< 



cC 



dist(z,S K )fc 2 / 3 ' 



Hence 



PII,1, X 

(z; k; t)(m p {z; k; t))' 1 -1+ "\ K ' 



in 



(3) 



2™ J go \(s)(s-z) 



-ds 



< 



cC 



dist(z, T, R )k 2 / 3 ' 



(10.69) 



(10.70) 



This is an extension of (5.35) to the case when z^0. 

For z = 0, from ( [10.63; ) and Ql0.64| ), R(0) = m< 3 )(0), and dist(0, E R ) > c x > 0. Note that A(s) is analytic 
in O except at s = — 1, and 



A(s) = -i2" 4 / 3 fc 1 / 3 [(s + 1) + ~(s + l) 2 + • • •], 



-1. 



By a residue calculation for ( |10.70| ), we have, as in (5.35) of |3j) 

^(.;M)W + jgga + 0(±). 



(10.71) 



(10.72) 



Thus using fl2.17|) and ( |2 . 1 8[ ) , from (gj), (|5.6|) , (jlO/Tj), ( |l0.8| ) and ( |l0.34D , we obtain Proposition |5T| (iii) of the 
case when < x < M. 



We now prove Proposition 5.2 when a; > 0. Since the choice of M was arbitrary in our calculations, for 
fixed x we choose M > large enough so that x < M. Let z E C \ £ be fixed. We first assume |z| < 1. 
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By modifying the contour T,( 3 \ if necessary, as in (10.50) and the following paragraph, we have z € Cl± and 
dist(z, E (3) ) > ci > 0. Thus from ( |10.70| ), \R(z) - I\ < Ck~ 1/3 with some constant C which depends on x. 
Thus from ( |5.6| ), ( |10.7 ), ( 10.34 ) and ( 10.64 ), we obtain the first limit of ( |5.13 ). Similar calculation applies to 
the case \z\ > 1, and we obtain the first limit of (5.14). The second limits of (5.13) and (5.14) follow from the 
first limits of ( 5.14 ) and ( |5.13| ) respectively by replacing z — > 1/z. Hence this extends the calculation (5.35) of 
[pi to the case when z is bounded away from the contour. 



Finally, we prove Proposition 5^ when x > 0. Set 

2 4 / 3 W 



a = 1 



fcl/3 



and 



fixed, 



x > fixed. 



2!/3fc2/3' 

In this case, —a € O. By a residue calculation again, for w not equal to 0, 

f 1 z2 4 / 3 1 

I 90 X(s)(s + a) S ~ (-l + a)fcV3 + aP^) 



1 

2?ri 



0( 



When it) = 0, we have the same order 0(k 1 / 3 ) by a similar calculation. On the other hand, since 



i-Vi-d) 



2t 
& 



fcV3 ^ 2 1 /3/ c 2/3 



using (|10.73|) we have dist(-a,S fl ) > Cfc _2/3 . Thus we obtain from ( |l0.70[ ) 

|JR(-a;fc;t)-7| < eC. 



Using A(— a) ~ —no, from (10.63) and (10.64), we have 

lim rrS 3 \~a;k]t) = m PII,2 (—iw,x) 

k — >oo 

since e is arbitrarily small. On the other hand, from the conditions on t and a, we have 

lim a k e- t{a - a ' 1) = el™ 3 " 23 ™. 

k — >oo 



Thus using (10.34), we obtain 

lim m (2) (-a; k;t) 

k — >oo 

lim mS ' \—a; k; t) 



m ' (—iw,x), 



m 1 ,2 (—iw, x) 



_g-|i£; 3 — 2xw 



(10.73) 



(10.74) 



(10.75) 



(10.76) 

(10.77) 

(10.78) 

(10.79) 

(10.80) 
(10.81) 



lim mS 2 \—a\ k; t) 

k — >oo 



m PII ' 2 (—iw, x) 



1 



-e 3 



8,„3 



w*+2xw 



-a e o 



(3) 



(10.82) 
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Now finally using (10.58), for each fixed w and x, we have 



lim rrv- 2 ' (— a; k; t) = m(— iw; x). 

k — >oo 



(10.83) 



From (10.7)-(10.10) and (10.79), this implies Proposition p.4| of the case when x > 0. This was a new compu- 
tation we had to do in the present paper in order to include the case when a — ► 1. 



10.2 When (2t)/k > 1. 

Throughout this subsection, we set 



It 



> 1. 



We need some definitions from 
(4.13) of §), 



Set < 9 C < it by sin 2 % 



(10.84) 



- . Define a probability measure on an arc (see 



dn(6) :=-«»(-) 



'--sin 2 (-)d0, 



-e„<e < 



(10.85) 



(10.86) 



(10.87) 



The measure dfi(0) is the equilibrium measure of a certain variational problem and the constant I is a related 
constant (see Section 4 of ||). For each \9\ < 9 C , the branch is chosen such that log(z — e l8 ) is analytic in 
C \ (—oo,—l] U {e 1 ^ : —tt < <fi < 9} and behaves like logz as z E K — > +oo. The basic properties of g(z) are 
summarized in Lemma 4.2 of [0. In general, the role of the ^-function in RHP analysis, first introduced in 
[ p79| ] and then generalized in p^7| , is to replace exponentially growing terms in the jump matrix by oscillating 



and define a constant (see (4.14) of [||) 

Z := — 7 + log 7 + 1. 
Now we introduce the so-called g-function (see (4.8) of ||), 

g(z;k;t) := log(z - e ,e )d/i(6»), z e C \ S U (-co, -1] 



or exponentially decaying terms. In [031, the authors introduced a g-function of a form similar to (10.87) to 



analyze an RHP associated to orthogonal polynomials on the real line. The above ^-function (10.87) introduced 
in H is an adaptation of their work to the circle case. When < 7 < 1, the related equilibrium measure is (see 
(4.12) of |) 



dn(0) = ^-(1 + 7Cos6»)d6>, 
2n 



-TT < 9 < TT, 



with the related constant I = 0, and hence (see (4.15) of 



;(*) 



logz 



SL 

2z ' 



-ttZ + iri, 



\z\ > l,z$ (-oo,-l) 
\z\ < 1. 



(10, 



(10.89) 



Since g(z) is explicit in this case, we did not introduce it of the form (10.87) in the previous subsection 
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Then up to ( |10.104| ), we follows the procedure in Section 6 of §]. Write S = C\ U C 2 where C 2 := {e ie : 
] c <9 < 9 C } and Ci := E \ C2. Define ( z; fc; t) by 



m (1) {z;k;t) := e% a3 Y(z; fc; t)e- ka{ -*' k ' t)a * e 
where 03 = (J Then solves (see (6.1) in ||) a new RHP 
m^(z; fc; t) is analytic in C\ E, 

(g— 2ka(z:k\t) t_ j\fc 

g2fcQ(z;fc;t) 
I (— l) fe e~ 2fe 5 — (z;fc;t) ' 

v° 1 

m^(z]k\t) = 7 + 0(i) as z -> 00, 

where 5(z;fc;i) is defined by (see Lemma 6.1 in ||) 

s + 1 



(10.90) 



(z;k;t) — m_ [z; fc; t) 



on C 2 , 



on Ci, 



(10.91) 



a(z; fc; t) := 



7 



■yj(a-e ie °)(s-e- ie °)da, £ := e iS " 



(10.92) 



(Notation: We use 5 here instead of a in || to avoid confusion with a in ( 10.140| ) below.) The branch is chosen 
such that ^/(s — e l9a )(s — e~ t8c ) is analytic in C \ C\ and behaves like s as s e M — > +00. 

Define m^ 2 ) (z; fc; i) as in ( |10.3| ). Thenm^ 2 ) solves a new RHP, normalized as z — > 00, with the jump matrices 
((6.2) of |) 



w( 2 )(z;fc;t) = 
v^{z-k;t) = 
Through the changes Y — » — > mfi\ we have 



I _ e -2ka(z;k;t) 
g2kai(z:k;t) q 
g— 2/cq — (z:k;t) ^ 

1 



on C 2 , 



on Ci. 



(10.93) 



Y u (z;k;t) = -e fcff( * ifcit) m^ («;*;*), 
y 21 (z;fc;i) = -(-l)V»<***>+ w m<?(z;M), 
y„(z;fc;t) = e fc »<**M?(*;M). 



|z|<l, 
|z| <1, 
|z| >1, 

w > 1. 



(10.94) 
(10.95) 
(10.96) 
(10.97) 



Set I]' 3 ) := C\ U Cinside U Coutside as in the Figure which divides C into four regions, f2*- , j = 1,- • • ,4. 
Again, there is a certain freedom choosing the shape of Ci ns id e and C ou t S ide- For example, Ci ns id e (resp., 
C ou tside) can be any smooth curve lying in fi 2 (resp. O3 ) connecting e l9c and e _ie<: ; the precise requirement 
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Figure 6: S^ 3 ^ and fi^when 7 > 1 



is given in g] (see also ( |10.100| )-( p.l02[ ) below). Define m^(z;k;t) by (see p.1151 of |3)) 



7.(3) = m( 2 ) 



i(3) = m( 2 ) 



7.(3) = m(2) 



2 p — 2fca(z;fc;i) 

,o 1 

1 0' 

a 2ka(z;k;t) j 



in fif, 



in £7 



(3) 
3 1 



)(3) n (3) 



in 



Then solves a RHP, normalized as z — > 00, with the jump matrix given by 



; (3) (z;fc;i) = < 




on Ci, 



on C, 



outside ) 



on Ci. 



From the properties of 5(2), it is proved in g] (between (6.4) and (6.4)) that 

e -ka-(z;k;t) _^ q ag _> ^ zeC"i, 



fca(js;/c;i) 



as fc — > 00, z G C, 



(10.98) 



(10.99) 



outside • 



(10.100) 

(10.101) 
(10.102) 



The choice of Ci ns id e and C ou t S ide is precisely for these properties. Here the convergence is uniform for any 
compact part of the each contour away from the end points e c and e~' l9c , but is not uniform on the whole 
contour. This gives rise the technical difficulty which will be overcome below using the idea of parametrix. 
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Formally v ^ — > i? 00 as — * oo where 

r 

V°°(z) 
v°°{z) 




On Cinside ^ Ccmiside 



on Ci. 



(10.103) 



Thus we expect that m^ 3 - 1 converges to m°°, the solution of the RHP mrj° = m^v 00 with to 00 — > I as z — > oo. 
The solution m°° is easily given by (see Lemma 6.2 in 0) 



(10.104) 



where /3(z) := ( z z _ e _ta c ) which is analytic C\ C\ and /3~+laszel^ +oo 



10.2.1 The case 2t > afc for some a > 1. 

The parametrix m p {z) is introduced in (6.25)-(6.31) in j^] , and has the following properties. In the neighborhood 
O of size e around the points e %e " and e~ l6,c , m p (z) is constructed using the Airy function in such a way that 



(m p (z)) + = (m p (z))^v^(z) for z S E (3) n O, and Hn^m 00 )" 



/I 



L°°(ao) 



O^" 1 ). In C \ O, we set 



m p (z) 



J (z). Then the ratio R(z) := m^(z)m p x (z) has no jump on E^nO, and has ajump 



W R - 



converging to I uniformly of order 0(k x ) on dO, and of order 0(e ck ) on E^ 3 ) n O c as fc — > oo. This implies 
that = 7 + O^ 1 ) for any z e C \ E p , E p := (E( 3 ) n O c ) U dO. Moreover following the arguments in 

Section 8 ofJb|, the error is uniform up to the boundary in each open regions in C \ E p . In particular, for 
z e n[ 3) U nf } , (see (6.34)-6.40) in @) 



m (3) (z) = (7 + 0(fc- 1 ))m oo (z). 



(10.105) 



Here the error is uniform for ak < 2t < bk for some < a < b. For the case (2t)/k — > oo, by shrinking the 
size of properly, we again obtain uniform error (see pi). Therefore for any a > 0, we obtain uniformity in 



(10.105) for ak < 2t. 



When z = 0, (3{0) = -ie l9 -/ 2 and g(0) = iri (Lemma 4.2 (vi) of ||), Also ( |10.98| ) says m( 3 )(0) = m^(0). 
Thus Proposition |t| (v) follows from ( |10.94| ) and ( |10.95| ), as in (6.4) of g. 

Now we consider Proposition 5.6 of the case where x < — 2 1 / 3 (a — l)fc 2 / 3 . For z = —a real, a)| = 1, so 



m°°(— a) is bounded. Hence from (10.98) and (10. 94), (10. 96), we have for a > 1, 

|7r fc (-a;fc;t)| < (7|e fcfl (-°"*^|. (10.106) 
Then we proceed as in ( 10.109| )-( |10.121 ) of the following subsubsection to obtain the proper estimate. 
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10.2.2 The case k + M 2- 1 / 3 k 1 / 3 < 2t < ak for some a > 1 and M > M . 

In this case, the points e ie " and e~ %0c are allowed to approach —1, but the rate is restricted: 

k \ 1/2 Jm 



H=2(l-^ >j^, (10-107) 



for k large. We now take the neighborhood O to be of size ^ — 1 around e c and e c . From (10.107), 
consists of two disjoint disks and their boundaries do not touch the real axis. We introduce the same parametrix 
m p as in the previous subsubsection. Then we have a similar result: there is Mq > such that for M > Mo, 

m(3)(z)= ( I+Q ( M |-l) )) mQO(z) - (10 - 108) 



for z e 0i 3) U «i 3) ; this is proved in (6.34)-(6.40) of 



When z = 0, as in the previous subsubsection, we obtain Proposition 5.1 (iv), as in 



Now we prove Proposition 5.6 when x < —M. As in the previous subsubsection, we have 

\TT k (-a;k;t)\ < c\e k9 ^ a *^\. (10.109) 

Now we need an estimate of g(— a; k; t); as we mentioned at the beginning of Section [l^, this is the second way 
in which we must extend || . Note that 

1 C e " 

Re(g(-a)) = - log(l + a 2 + 2a cos 9)dfi(9) 

2J -0* (10.110) 

= log2 + - log -+/(s), 
2 7 

where 

I(s) = - [ log(s 2 - x 2 )y/l-x 2 dx, s := 1l a) > 1- (10.111) 

The inequality s > 1 follows from the arithmetic-geometric mean inequality and the assumption 7 > 1. A 
residue calculation gives us 

I'{V) = ~I "i^Vl - x 2 dx = 2y- 2^/i/ - 1, y>l. 

t J-i y - x 

Integrating from 1 to s > 1, we have 

I(s) = s 2 -1-2 J Vy 2 - ldy + 1(1). 
The constant 1(1) can be evaluated (cf. Lemma 4.3 (ii)-(a) in ||3) ) : 

7(1) = - I 2 log(sin 2 0) sin 2 0d(9 = ^ - log 2. (10.114) 

7T 2 



(10.112) 



(10.113) 
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Thus we have 



Re(ff(-a)) = -- + - log - + s 2 - sV ' s 2 - 1 + log(s + y/s 2 - 1). 



Assume < a < 1. We change the variables 7, a into s, £ where s is defined in (10.111) and 

1/2 



> 1. 



Then 



F(0 := ff(-a) - |a = -i - log^ - ^ 2 + 2s£ - s 2 - sVs 2 -l + log(s + - 1). 



(10.115) 



(10.116) 



(10.117) 



Differentiating with respect to £, we find 

F'(0 = -i-£ + 2s. (10.118) 

Thus the maximum of F occurs at £_ = s + sjs 2 - 1. But F(s + Vs 2 - 1) = 0, hence F(£) < 0. Thus we obtain 

\e- ta Ti k {-a- k)\ < Cf e kB»(fl(-o:*;*)-*o) < C) < a < 1. (10.119) 

For a > 1, note that 

Re(>(-a)) = loga + Re(.g(-a" 1 )). (10.120) 



Thus using (10.119), we have 

|e-* a_1 (-a)- fe 7rfe(-Q!;fe)| < c r e fcRe ( 9( - a;fc; * ) -* a " I - loga ) = c , e fcRx! ( s( - Q ~ 1;fc;t) -2 Q ~ 1 ) < c. (10.121) 

10.2.3 Case k < 2t < k + M2~ 1 / 3 fc 1 / 3 for some M > 0. 

First we introduce m PI1 - 3 as in (2.22)-(2.28) of @. Set 



(10.122) 



which is analytic in C\ [— y y -jp] and behaves like ^z 3 + xz + ^ + 0(z 3 ) =: 6pjj(z) + 0(z 1 ) as z — > +00. 
Let £ p// < 3 := U^E^ 7,3 as shown in Fi gure f[| The angles of the rays with the real line are between and 



7r/3. Recall that m(z;x) solves (2.15) the RHP for PII equation. Define m PII,3 (z;x) by 



,PII,3 _ 



m(z; x)e 



Ka -8pii)<?3 



P/J,3 P//,3 



= m(z; x) 



m(z; x) I j ' 



1 

1 (i ' 



(10.123) 



-,P//,3 , , r ,PII,3\ 



Z G (^2 U 
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PII.3 PII.3 




Then mS 3 ' solves the RHP (see (2.25) in H) normalized at oo with the jump matrix 



v (3 \z;k;t) 



1 ( ) 

P I I 



1 

_ e -«fl* 
1 





^PII,3 ^P//,3 

on S 1 ,E 2 



^PII,3 ^PII,3 

on E 3 ,S 4 



on £ 



PII.3 



1 



Also we have 



where m PII,3 (z; x) = I H — 



Now as before, set x by 



m 1 {x)=m 1 {x)~{ — )a 3 . 



+ 0(z 2 ) as z — > oo. 



2f 



l 



2 l/3 fc 2/3- 

Hence we have — Af < a; < in this subsubsection. We now proceed as in case (iii) of Section 6 in 
the parametrix 



m p (z;k;t) = m P7/ < 3 (A(z), a;) 
m p (z; k;t) = I 



in e>\£( 3 ), 
in O c \£( 3 \ 



(10.124) 



(10.125) 



(10.126) 
|. Define 

(10.127) 



where X(z) is defined in (10.62) and O is a small neighborhood of size e > around z = — 1 (sec H case 
(iii) of Section 6 for details). As in Subsubsection 10.1.3 , the ratio R(z; k;t) :— m^nip 1 satisfies a new RHP, 



normalized at oo, with jump matrix vr satsifying the estimate (10.66) where m PII ' 2 (x) is replaced by m J ^ 11 '' > (x). 
Hence we have 



PII.3, 



m {3 \z;k;t)(m p (z;k;t)) - 7 + — ^ 



PII,3 



(x) 



1 



27T« J d0 \(s)(s~z) 



-ds 



< 



cC 



dist(z, E fl )fc 2 / 3 ' 



(10.128) 



which is hidden in the derivation of (6.19) of 
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Then as in (10.72), we have 



which is a (direct) extension of (6.19) in §]. Now from fll0.98j ) and ( |10.125[ ) (see (6.19) of (|), 

m (2) (0;fc;<) = /+ u '— - + 



k l/3 



k 2/3, 



(10.129) 



(10.130) 



Hence using e kl = 1 - 25/ f^ 1/3 + 0(k 2/3 ) and .g(0) = ni, ( |l0.94| ) and ( |10.95|) yield Proposition |1] (iii) in the 
case — M < x < 0, as in ||. 



For the proof of Proposition bJ2, note that, as before, for each fixed z S C \ S, we can use the freedom of 
the shape of X^ 3 ' (and £#) so that z G 17 ^ U fl^ , and dist(z, S^) > ci > 0. Thus we obtain 



lim m (2) (z; fc; t) = /, zeC\S fixed. 



(10.131) 



From (10.94) and (10.96), we have 



lim e- k ^ z ' k ^TT k (z) = 0, \z\<l, 



lim e- fe9 ^ fc; *V fc (z) = 1, 

— >oo 



z > 1. 



(10.132) 
(10.133) 



This is an extension of the calculation of |3| where z = is given. Now in order to prove Proposition 5.2 



wo 



need further analysis which is an extension in the second category as mentioned as the beginning of Section 10 



Since 7 = (2t)/k, for the proof of Proposition 5.2, it is enough to show that 



lim (_i)fc e M9(z;fc;t)+^] = 1; \ z \ < ^ 

k — >oo 

lim e *fa(*;*;t)+3*- l -iog*] = 1 ui > L 



(10.134) 
(10.135) 



But the proof of Lemma 4.3 (ii) of M says that for \z\ > 1, z ^ (—00, — 1), 



.9(2) = T^log^- + z ^ + ^+4 



7 , 7 



(s - e iB °)(s - e- ie -)ds +#-(!), 



(10.136) 



where the integral is taken over a curve from l + o to z lying in {z € C : |z| > 1,2 ^ (— 00,— 1)}. Here 
a/(s — e l9c )(s — e~ l9c ) is analytic in C \ C\ and behaves like s as s G R — > +00, and logz is analytic in 
C \ (—oo,0] and is real for z]R + . Calculations in the same proof, together with Lemma 4.2 (viii) of H, give us 
9-{l) = ~\- 5l°g7- Also using sin 2 % = A for |s| > 1, s (-00, -1), 



(s - e ie -){s - e- ie °) = (s+ 1) - 



2s 
s+1 



(7-l) + 0((7-l) 2 )- 



(10.137) 



Thus expanding in 7 — 1, we have 



g(z) + jz- 1 - logz = 0(( 7 - I) 2 ) = 0(fc- 4 / 3 ), 



(10.138) 
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which implies (10.135). Similar calculations using for \z\ < 1, z ^ (—1,0] 

s + 1 



<?(z) = ilogz- J(z + z + | + J 



y(s - e l9 -)(s - e-' l9 -)ds + g+{l), 



(10.139) 



and g + (l) = -\ - \ log z + iri yield ( [UUgJ ). 



For the proof of Proposition 5.4 when x < 0, set 



a = 1 



2 4/3 



fcV3 



(10.140) 



In this case, we need a new argument as a — > 1. When to and # are fixed, again as in ( 10.771) , we have 
limfc^oo R(—a; k; t) — I, which implies that 

lim m (3) (~a;k;t) = m PII ' 3 (-iw,x). (10.141) 

k — >oo 

From (6.8) of we have 

A - TQ/n 

(10.142) 



Jim fca(— a; fc; t) = i— ((— iw) 2 + ^) 3 ^ 2 = ig PH 



which from ( |10.98|) and (|10.123| ) implies 

lim m^(-a;M) = m(-no, a ;)e < k P "(- <to )- fl ''"(- <,0 >> \ (10.143) 

— >oo 

Now we compute the large k limit of kg{—a; k; t) — ta when w > 0, and of kg{— a; k; t) — ta _1 — log a when 

r sin 8 
^ a+cos t 



w < 0. For — 7r < 9 < 7r, lim e jo arg(— Q + «e — e ie ) = 7r + tan 1 ( - "" Q g g ) where — 7r < tan 1 <f) < it. Since 
tan-^ ^^g ) is odd in 9, we have from (|l0.115| ), 

limj(-a + ie) = lim Re g(— a + ie) + ni 

= + X lo S - + s 2 - s Vs 2 - 1 + log(s + Vs 2 - 1) + 
2 z 7 

where s = ^~|yi=~^ > 1. Under the stated conditions on 7 and a, as k — > 00, 



(10.144) 



- 



V7(l + a) = / ^ 2 

2VS U 1/3 2 4 / 3 7 fc 2 / 3 ' k 

2^w 



1 2.3 Q( x 



7 



1 



2xu> , 1 s 



(10.145) 



(10.146) 



fcl/3 2V3fc2/3 fc V /c-4/3- 

Note that lim^o g{— en + ie) — \im £ io g(—a — ie) is 2iri for a > 1, and is for < a < 1. Therefore we obtain 



lim (_i)* e M-°0-ta = e s 



i ffi 3_„_l / , lt 2_.^/2 



-f(^-f) J 



lim (— l) fc e fc 9( _a ) _tQ 1 -k\o S a 

k — >oo 

Also being careful of branch, we have 

i(g PI1 \—iw) — Opii(-iw)) = 

i{g PII {-iw)-e PII {-iw)) = 



4 3 4 f 2 x \ 



3/2 



— 10 — xw + — 1 w — — 



3/2 



2 



w > 0, 
10 < 0. 

to > 0, 
w < 0. 



(10.147) 
(10.148) 

(10.149) 
(10.150) 



Since lim^^ e H = 1, using ( |10.94| )-( prl97| ) and ( ]10.143| ), this implies (|5.21|)-(|5~24|) when x < 0. 
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10.3 Proof of Proposition |5.8 



The analysis in this section is new, and is needed for the proof of Proposition 5.8 . Let a > 1 be fixed. Set 



Q,f a 2 _ 1)1/2 

-==, x e R \ {0} is fixed. 



* 2 + 1 (a 2 + l) 3 / 2 Vfc' 



(10.151) 



We are interested in the asymptotics of e Qt (— a) l irk(— a 1 ;t). Since a 2 ° 1 < 1 and a; is fixed, we are in the 



case of Subsubsection [10. 1.1] and/or 10.1.2. We define and as in (10.1) and (10.3). Recall that we 
have certain freedom in the choice of Y,( 3 >. We choose a contour passing through the saddle points of 



f(z) :=-( z - z -^+log(-z), 



(10.152) 



the exponent of the (12) entry of divided by k. The saddle points (see ( 10.31 ) and the following discussion) 
are — p v and —p^\ where 



1 - y/1 - (2i/fc) 2 1 (a 2 + l)!/ 2 x 1 l 



(10.153) 



Take <5 > and e > small such that E c := {— p c + is : —k s 1 ' 2 < s < e/c 5 1 / 2 } lies inside the open unit disc 
for all fc > 1. Define (see Figure |) S< 3 ' := E^ U by E,^ } := E c U E r and E^ t : = {r"V* : re** € E^ } }, 




Figure 8: E< 3 ) and fi< 3 ) 



where E r := {| - p c + iek 5 - l ' 2 \e i6 : \9\ < 9 Q }, -p c + iek 5 - 1 ' 2 = | - p c + iek 5 ' 1 / 2 ^ . Let Vtf ] be as in Figure 



| and define rrS^ as in ( 10.34 ). As in ( 10.46 )-( 10.48 ), the quantities we are interested in are 

e -«*(_ a )fe 7rfe (_ a -i ; t) = -o/'e-^ a - cr ^rt^{-a- 1 \k]t), x < 0, (10.154) 

e - a *(-a)*7rfc(-o -1 ;t) = -a k e- t ( a-0 " I )mS ) (-a _1 ;*;t)+mg ) (-a- 1 ;fe;t), J > 0. (10.155) 

For the estimates of u/ 3 ) := — I, note that for any fixed < p < 1, Ref(pe td ) — F(pe %9 ; ^) (recall 
10.11)) is increasing in < 8 < n, and is decreasing in tt < 9 < 2ir, hence ||e fe ^ z )||£oo(E r ) — e fcRc -^~ Pc+efc ' 5 1/2 ) _ 
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But we have 



f(-Pc +ia) = ^(a-a x ) - log a - ^ f 



1 /a 2 (a 2 - 1) 



a 2 + l 



+ 0(k-l +2S ), \a\<ek s - x /' 2 . (10.156) 



Thus, 



(10.157) 



and also using 



-k p t(a-a-i) = p *[l°e«-^T+o(^)] 



log a 



a 2 -l 



> 0, for a > 1, 



(10.158) 



we have 

l|e fc/(z) l|L~ (Sr ) <Ce- cfe . (10.159) 

On the other hand, one can directly check that Ref(—p c + a) has its maximum at a = for — e/c 5 " 1 / 2 < a < 
efc 5 - 1 /2 ; ae nce ||e fc/ ( 2 )|| L = (Ec) = e kRcf( -~ p "l Again ( |10.156j ) and Q10.158| ) yield 

||e fc/(z) || L ~ (Ec) <e- ck . (10.160) 

Similarly, we have ||e _fe ^^^|| roo , v (3) < e~ cfc . Now calculations as in Subsubsection 10.1.2 give us the result 

( |10.43D . Hence using ( |10.158| ) and ( |10.160| ) and noting dist(-or\ £^) = ^rjyr^ ■ t^, we have 



,(3) 



(10.161) 



ake -t( a - a -i) m (3) 



m\ 2 ' (-a ;k;t) 



Q fe e -t(a-a- 1 ) 



(- fl )fc e *('-" 1} 



2?ri 



0(Vfce 



ce ™ j. (10.162) 

To evaluate the integral asymptotically, first we change the variable by s = — p c — ^(a 2 —!) 1 / 2 ' vT' 



• -H=. Then 



from (10.156), the numerator of the integrand becomes 



a -fc e i(«-a-%-A(yW)+0(fc-2 +25 )_ 



Hence setting .4 



._ taja 2 -!) 1 ' 2 
( Q 2 + 1 )l/2 : 



(10.163) 



1 

2ni 



(-s) fc e 



k t(s-s x ) 



-ds 



1 



27ri J-Ak s y + ix 



dy(l + 0{k- 2+2A )) + 0(e~ ck ). (10.164) 



Thus from ( |10.154[ ) and ( |10.155[ ), we obtain 

\k 



lim e- at {- a )«-K k {-a- x ;t) = — 

k— >oo Z7TI 



1 f°° e -hiv +x ) 



lim e" Qt (-a) fe 7r fc (-a- 1 ;i) 



l 

2?ri 



e -±(j, 3 +s 3 ) 



J/ + M 



dy + l, 



x < 0, 
x > 0. 



The function h(x) :— ^ e \ + i X — dy is smooth in x > and x < 0. The derivative is h'(x) — "^^ e 
As x — > ±oo, /i(x) — > 0. Therefore we see that 



(10.165) 
(10.166) 



1 f 00 e ^^ y2+x2 " 



2iri 



y + ix 



-dy 



1 fl — i-U 2 

-7= c 

V27T J oo 



J x e 2V dy, x > 0. 



(10.167) 



Thus we proved Proposition 5.8 
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